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ABSTRACT 

Acoustic  waves  of  finite  extent  were  propagated  into  a  dc  glow 
discharge.   It  was  found  that  compression  phases  of  these  waves  reduced 
light  emitted  by  the  discharge,  while  rarefaction  phases  increased  this 
light.   In  the  compression  phase  of  the  acoustic  wave  these  changes  were 
attributed  to  increased  losses  of  electron  energy  through  an  increase  in 
the  collision  frequency.   Electron  energy  was  increased  in  the  rarefaction 
phase  of  the  acoustic  wave  as  collision  frequency  was  reduced. 

Fast  and  slow  waves  of  stratification  were  excited  by  the  acoustic 
waves.   These  waves  of  stratification  had  opposite  signs  of  light  in- 
tensity variation  when  excited  by  compression  and  rarefaction  phases  of 
the  acoustic  wave. 
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CHAPTER  I 
INTRODUCTION 

Electric  and  magnetic  fields  have  been  widely  used  for  investigation 
of  the  properties  of  glow  discharges  because  of  their  strong  interaction 
with  charged  particles.   In  recent  years,  the  interaction  of  acoustic 
waves  with  the  discharge  plasma  has  received  some  attention.   These  ex- 
periments have  been  limited  in  scope  and  little  theoretical  work  has 
accompanied  them.   This  work  is  an  experimental  study  of  the  effects  of 
an  acoustic  wave  group  of  finite  extent  upon  a  neon  glow  discharge  40  cm 
long,  in  a  2.5  cm  diameter  tube  in  the  pressure  range  0.7  to  4.0  torr, 
and  operated  at  currents  up  to  35  ma. 

In  the  discharge  of  electricity  through  low  pressure  noble  gases, 
spontaneous  oscillations  of  the  discharge  current  are  nearly  always  ob- 
served.  Frequencies  of  these  oscillations  range  from  0.5  to  50  kHz. 
Photomult iplier  signals  indicate  that  the  light  intensity  of  the  positive 
column  also  oscillates  with  this  frequency.   These  oscillations  are 
caused  by  the  anode  to  cathode  motion  of  light  areas  in  the  discharge 
which  are  known  as  moving  striations. 

Studies  of  moving  striations  have  made  great  progress  in  recent 
years  by  investigating  transients  introduced  into  a  homogeneous  discharge 
through  impulsive  modulation  of  the  current,  or  through  application  of  a 
high  voltage  pulse  to  a  ring  encircling  part  of  the  positive  column. 
These  transients  have  the  form  of  a  wave  group  or  wave  packet  which  prop- 
agates toward  the  anode  with  a  group  velocity  u,  while  the  component 
waves  have  phase  velocity  v  directed  toward  the  cathode.   It  is  therefore 
a  backward  wave.   Individual  waves  have  the  characteristics  of  moving 
striations,  and  the  group  has  been  termed  a  "wave  of  stratification". 


Figure  1-1  shows  schematically  the  traces  observed  by  a  photomulti- 
plier  at  two  points  in  the  column  (1-1 (a),  (b)).   From  these  the  propaga- 
tion in  the  direction  of  the  anode  can  be  seen.   Part  (c)  of  this  figure 
shows  development  of  the  wave  of  stratification  in  time  and  space,  with 
positive  X  chosen  in  the  direction  of  the  electric  field.   In  this  dia- 
gram, the  dark  areas  correspond  to  the  moving  striations  which  are  present 
within  a  wave  packet.   Striation  motion  toward  the  cathode  is  evident, 
and  the  striation  amplitude  is  depicted  by  the  width  of  the  line. 

Under  conditions  in  which  this  wave  receives  large  amplification 
during  propagation  along  the  column,  the  oscillations  of  current  which 
occur  as  it  arrives  at  the  anode  cause  excitation  of  additional  waves  of 
stratification.   Further  increase  of  the  amplification  causes  the  striations 
to  be  self-excited,  with  striations  of  largest  amplitude  appearing  at  the 
anode  end  of  the  discharge. 

Pekarek  first  reported  these  waves  in  1954,  and  his  ensuing  theoret- 
ical and  experimental  works  have  established  several  important  character- 
istics :  [  41  I 

1)  the  waves  are  ionization  phenomena,  dependent  upon  the 
change  of  ionization  rate  with  electron  temperature; 

2)  a  single  slow  and  two  fast  waves  exist,  with  the  slow  wave 
dependent  upon  stepwise  ionization,  the  fast  waves  on  direct  ionization; 

3)  the  slower  of  the  fast  waves  is  associated  with  motion  of 
an  atomic  ion,  the  faster  with  the  motion  of  a  molecular  ion; 

4)  any  localized  perturbation  of  ionization  rate  of  sufficient 
amplitude  will  cause  the  wave  to  originate;  perturbations  of  opposite 
sign  produce  waves  with  opposite  courses  of  light  intensity; 

5)  the  formation  of  space  charges  due  to  higher  mobility  of 
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Figure   1-1.       (a)    Time   modulation   of    light    intensity    in   anode 
region  by   wave   of   stratification. 

(b)  Time  modulation   of    light    intensity   nearer 
cathode. 

(c)  Time-space   development    of    light    intensity 
by  wave   of    stratification. 


the  electrons,  and  the  consequent  change  of  electron  temperature,  cause 
the  wave  to  propagate; 

6)  striations  with  small  separation  (or  wavelength,  J^   )  are 
damped  by  longitudinal  ambipolar  diffusion,  those  with  long  wavelength 
by  relaxation  of  electron  temperature. 

Pekarek  has  shown  that  any  spatial  perturbation  of  the  steady-state 
ion  density  causes  a  localized  space-charge  region  to  be  formed.   53 
On  this  basis  the  formation  and  motion  of  striations  can  be  explained  by 
the  following  chain  of  events: 

a)  ion  density  is  increased; 

b)  a  space  charge  region  is  formed; 

c)  electric  field  is  increased  on  the  anode  side  of  this  region, 
decreased  on  the  cathode  side; 

d)  electron  temperature  changes  with  the  electric  field,  causing 

e)  increase  of  light  toward  cathode,  decrease  toward  anode; 

f)  ions  move  toward  the  cathode  under  the  influence  of  the  main 
electric  field. 

Propagation  of  the  wave  toward  the  anode  can  be  explained  by  considering 
the  decrease  of  electric  field  on  the  anode  side  of  this  disturbance. 
In  this  case: 

g)  electron  temperature  and  ionization  rate  are  decreased, 
causing 

h)  a  reduction  of  light  intensity  and  a  negatively  charged 
region; 

i)  this  region  of  negative  charge  enhances  the  electric  field 
on  its  anode  side,  causing 

j)  increase  of  ionization  rate  and  ion  density. 
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Step  (j)  was  the  initial  condition  of  the  disturbance.   Since  the 
first  ion  density  perturbation  caused  a  second  similar  region  to  form 
on  the  anode  side,  repetition  of  this  cycle  describes  propagation  of  the 
wave  in  the  direction  of  the  anode. 

Properties  of  the  wave  can  be  related  to  microphysical  processes  in 
the  discharge  through  the  relaxation  time  which  is  defined  as 

r  =  -i-  •  — - —  en 

1  o  2  u  +   v  K    ' 

This   relaxation  time   is   characteristic   of  the   slowest   process    in  the 
chain  of  events   which  form  and   propagate  the  wave. 

In  the  case  of  the   fast   wave,    it    is    the  time   required   for  the   excess 
space  charge  to   vanish  through  diffusion   losses,   while    in  the  slow  wave 
it    is   the   sum  of  this  diffusion  time   and   the  time   required  to    ionize  the 
excess  of   ions    produced  by  the    increase  of  electron   temperature. 

Dispersion  of  this   wave   in  the    inert   gases    is    such  that  the  higher 
frequency  components   of  the   packet   generally   appear  earliest   at    a  given 
position.      In  time-space   photographs,    this   type  of  dispersion  causes   the 
striation  paths   to   be  curved   convex  with  respect   to   origin. 

Amplification  of  the  wave   is   dependent   upon  the   change  of    ionization 
rate  with  temperature,    and  upon  the  spatial   phase  difference  between  the 
course  of   electric   field  variation  and  that   of  the  electron  temperature. 

M 

Thus   the   essential   properties   of    the  wave  of   stratification  are 
well   understood.      Detailed  mathematical    formulation  of  the  wave  of   strati- 
fication has    been   carried   out   by   Pekarek   and   various   co-workers      41   -    60    , 
and  the  theory   for  the   case  of   zero  Debye   length  has   been   found  to  des- 
cribe this    case  reasonably  well.      However,    all   perturbation  work  has 
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been  accomplished  electrically  with  the  exception  of  some  work  by  Sicha 
et .  al.,  who  used  square-wave  modulation  of  the  power  input  to  re-entrant 
RF  cavities  to  excite  the  wave  of  stratification,  thus  directly  per- 
turbing the  ion  density  locally.  J  7 3  j  Other  areas  of  the  column  are 
affected,  however,  as  standing  striations  form  between  the  point  of 
application  of  the  RF  discharge  and  the  anode.   These  are  similar  to 
those  found  at  the  cathode  end  of  the  positive  column.   Both  the  high 
voltage  pulses  and  pulse  modulation  of  the  current  also  cause  disturbance 
of  the  discharge  over  large  areas. 

There  has  been  no  reported  attempt  to  excite  waves  of  stratifica- 
tion by  magnetic  fields,  by  the  use  of  radiation,  or  by  acoustic  waves. 
In  the  acoustic  work  which  has  been  performed  in  dc  discharges,  Suber- 
tova  has  reported  creating  standing  striations  by  the  use  of  acoustic 
waves,  80   Carretta  and  Moore  observed  acoustic  waves  which  they 
attributed  to  moving  striations   9  ,  although  these  were  later  shown 
to  originate  in  the  cathode  region  of  the  discharge,  and  attributed  to 
changes  in  gas  heating  as  discharge  current  was  changed.   83   Saxton 
investigated  the  influence  of  acoustic  waves  on  the  positive  column  and 
found  that  the  amplitude  of  the  electron  density  variation  was  propor- 
tional to  the  neutral  gas  density  variation.   11 ,  12   Wojaczek  observed 
that  a  pulse  discharge  at  one  end  of  the  main  discharge  caused  acoustic 
waves  which  modulated  the  intensity  of  the  light  emitted  from  the  posi- 
tive column.   87 

With  the  exception  of  a  brief  explanation  by  Saxton,  there  has  been 
no  analysis  of  the  acoustic  interactions  by  the  authors  of  these  works. 
Several  other  experiments  involving  acoustic  waves  have  been  performed, 
but  none  have  involved  dc  discharges.   Ingard  has  theoretically  analyzed 
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the   possibilities   of  producing   acoustic   waves   from  interaction  with 
charged    particles,    but    has    given    little   consideration   to   the   inverse 
problem.      27 

The  way   is   clearly  open  for   an  extensive   investigation  of  the  effect 
of   acoustic   waves   on  the   positive   column  of   a  dc  discharge. 

In  this   work,    a  pulse  discharge   is   used  to  excite   acoustic  waves 
of   finite  extent.      Amplitude  of  these  waves    is  measured   by   a  thermionic 
gas-filled  diode,    and    interaction   in  the  positive   column  of  the  dis- 
charge   is  observed   with  a   photomult iplier  tube  and   displayed  on  time- 
space   photographs.       14,7  6 

A  theoretical   explanation  of   the  interaction  of  acoustic  waves    in 
the   positive  column   is   given. 
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CHAPTER  II 
EXPERIMENTAL  APPARATUS 

Figure  2-1  shows  the  arrangement  of  equipment  used  in  this  experi- 
ment.  Acoustic  waves  were  excited  by  a  discharge  operated  between 
electrodes  El  and  E2 ,  and  observed  in  the  dc  discharge  between  electrodes 
E3  and  Ek.      The  exciting  discharges  were  produced  by  high  voltage  pulses 
in  most  cases.   Square  wave  modulation  of  the  discharge  current  was  used 
to  compare  the  effect  of  propagating  a  compression  or  a  rarefaction 
phase  as  the  leading  component  of  the  wave.   Hollow  cylinders  with  dia- 
meter 18  mm  were  employed  to  take  advantage  of  the  higher  efficiency  and 
higher  current  densities  which  characterize  the  hollow  cathode  effect. 
62,  63,  78,  79   Cross-sectional  area  is  small,  giving  minimum  inter- 
ference with  the  propagation  of  acoustic  waves.   Electrode  El  was  always 
used  as  the  cathode  of  the  exciting  discharge;  polarity  of  the  electrodes 
in  the  observing  discharge  could  be  reversed  easily  to  study  effect  of 
direction  of  propagation  in  the  discharge. 

Power  for  the  observing  discharge  was  provided  by  a  0-1  kv  supply, 
the  working  load  was  a  pentode  with  indirectly  heated  cathode  (RCA  807). 
This  was  chosen  to  give  ripple-free  operation  of  discharge  current. 
Variable  resistors  were  used  in  the  cathode  circuit,  allowing  the  grid 
to  be  biased  negatively  by  current  flow  and  the  operating  current  to  be 
selected  by  this  simple  change  of  resistance.   A  pentode  was  selected  as 
the  discharge  load  since  this  has  been  shown  to  increase  the  range  of 
current  in  which  the  discharge  is  striation-free.  \^\     Figure  2-2  shows 
the  current  and  pressure  boundary  of  striation-free  operation.   A  curve 
calculated  from  the  empirical  relations  in   69   is  also  shown.   For  a 
given  pressure,  both  the  upper  and  lower  current  boundaries  of  this  region 
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were  found  to  be  larger  than  the  values  given  by  the  empirical  relations. 
Formulae  for  these  relations  are  given  in  Appendix  II. 

For  observing  the  interaction  in  the  discharge,  the  time-space  dis- 
play method  described  by  Stirand,  et .  al. ,  was  employed  (Figure  2-1).  j  7  6  1 
In  this  technique,  the  oscilloscope  horizontal  sweep  is  triggered  by  the 
exciting  discharge,  vertical  position  of  the  trace  corresponds  to  posi- 
tion of  the  photomultiplier  along  the  column,  and  the  intensity  of  the 
beam  is  modulated  by  the  photomultiplier  signal.   With  the  carriage  at  a 
given  position,  variation  of  intensity  of  the  trace  indicates  a  time- 
varying  phenomenon.   When  the  carriage  is  moved  to  a  new  position,  the 
vertical  position  of  the  trace  is  changed.   These  variations  occur  at 
different  times  for  phenomena  which  propagate  with  finite  velocities.   A 
photograph  taken  with  the  shutter  open  while  the  carriage  is  moved  slowly 
along  the  column,  records  the  development  in  time  and  space  of  time- 
varying  phenomena.   Figure  2-3  is  a  time-space  photograph  showing  the 
propagation  of  an  acoustic  wave  from  anode  to  cathode  in  the  discharge 
tube.   The  anode  is  at  the  upper  limit  of  the  trace,  the  cathode  at  the 
lower  limit  and  time  increases  to  the  right.   Light  intensity  variations 
occur  earlier  at  the  anode  end  of  the  discharge  because  of  the  anode  to 
cathode  propagation.   Photomultiplier  output  is  at  negative  potential, 
causing  photographs   taken  in  this  manner  to  show  light  intensity  varia- 
tions in  an  inverted  manner,  i.e.,  increases  in  light  intensity  in  the 
discharge  correspond  to  darker  areas  in  the  photograph.   The  unmodulated 
level  of  illumination  of  the  oscilloscope  screen  was  chosen  such  that 
very  small  variational  signals  could  be  observed.   Frequency  response  of 
the  intensity  modulation  circuit  caused  distortion  of  light  intensity 
changes  longer  than  about  five  milliseconds. 
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Density  variations    in  the  acoustic  waves  were  measured  by   recording 
changes    in  the  conductivity  of   a  gas-filled  thermionic  diode.      This 
technique,    first    noted   by  Dayton,    et .    al.       14 1,    is    based   upon  mobility- 
limited  operation  of   the  diode  according  to  the  equation 

i   =  KV2.  (2) 

Here  K  involves  the  mobility  which  is  a  function  of  neutral  particle 
density.   This  device  is  calibrated  by  recording  current  as  a  function 
of  pressure  with  voltage  fixed.   In  this  case,  the  density  changes  are 
directly  related  to  pressure  changes  through  the  perfect  gas  law,  as 
temperature  is  also  held  constant.   The  functional  dependence  of  current 
upon  pressure  can  therefore  be  expressed  as 

i  =  -^-   ,  (3) 


_/>m   ' 


where  the  value  of  m  is  obtained  from  the  curves  of  figure  2-4.   Fraction- 
al density  variations  result  in  fractional  current  changes,  which  may  be 
read  as  voltage  variations  across  a  resistor  R: 

P  miQ     m.R.iQ 

From  the  data    in   figure   2-4,    for   a  voltage  of   20   volts   m   is   0.286, 
giving 

-*£  =    3.47      «-   ,  (5) 

-r  o 

with   R  =    1000   ohms,    and    i      in  ma. 

In  the   acoustic  waves,    compressions    and  rarefactions   are  adiabatic, 
with  pressure   and  density   related   through 


P 


yy 


=   const.  (6) 
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Figure  2-4.   Diode  current  (multiplied  by  2)  as  a 

function  of  pressure  for  constant  plate 
vo 1 tage. 
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Thus 

_ip_  =  A?  .  y  (7) 

p        " 

and   the   sensitivity  of   the  device   at   1.0   torr,   with  20  volts   from  anode 

to   cathode   is 

S=— ^—  =0.24  ^—  .  (8) 

dp  x<<bar 

Dayton's  method  gives   a  sensitivity  of  0.79  mV^ubar   in  nitrogen  at    1.0 
torr.      This   can  be   attributed   to   the   larger  total   collision  cross-section 
of  nitrogen.      His   calculation  may  be   in  error,   however,    as   he  seems   to 
have   assumed  an   isothermal   relation  between  density  and  pressure   instead 
of   an   adiabatic   relation. 

Experiments   showed  that   operation  of  the  neon-filled  diode   in  the 
experimental   pressure  range  was   described   by  the  equation 

i   =    KV1*6    .  (9) 

This    indicates   that   the  diode   is   nearly   space-charge   limited   rather 
than  mobility   limited.      The  curve  of  figure   2-4   indicates   that   the  varia- 
tion of  conductivity  with  pressure    is   sufficient   to   allow  its   use  as   a 
detector  of  density  variation. 

Acoustic   wave  amplitudes    are  expressed   in  fractional  variations  of 
neutral   gas   density,   as    it    is   these  variations  which  will   affect   the 
properties   of  the  plasma.      All   measurements  of   acoustic  wave   amplitudes 
were  made  with  the  observing  discharge   extinguished. 

Figure   2-5    shows   the  diode   signals   for  acoustic  waves   produced   by 
the  high  voltage  pulse  discharge   and  by  square  wave  modulation  of  the 
exciting  discharge  current   by  the  pentode.     Density   increases   are   indi- 
cated by   a  negative  change   in  voltage.      In  this    figure   it    is   evident 
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5   msec/division 

Figure  2-5.       (a)   Acoustic  wave  produced   by  2   kv  on  plate  of 
high  voltage  pulser.      Density  variation   is  0.27   per  cent 
(b)   Acoustic  wave  produced  by  square  wave  modulation  of 
discharge  current.      Peak  density  variation   is  0.094  per 
cent,  difference   in  static   levels    is   0.04  per  cent. 
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that  the  high-voltage  pulses   excite  waves  of  much   larger   amplitude  and 
shorter  rise  time.      Maximum  density  variation   is  0.27   per  cent   for   the 
larger  amplitude  acoustic  wave,    and   0.09  per  cent   for  the  smaller  am- 
plitude wave.      A  static  density  variation  of  0.04  per  cent    is   also 
created  by   square  wave  modulation  of  the  discharge  current. 

Translation  of   the  gas   along  the  discharge  tube  axis   by  a  sinusoidal 
acoustic  wave   is  given  by 

X 


i'l 


C.S0«SinCUt  dt 
o 


CS 
=  ~^f-   (1   -  cosuit)    ,  (10) 

where 

C        =        sonic  velocity      (456  m/sec  @  24°C) 

S        =        condensation    (df/p). 

For  the  high  voltage  pulser,    a  frequency  of    1   kHz   approximates   the  com- 
pression phase,    and  the  maximum  translation   is   1.23  mm  when  the  density 
variation   is   0.27   per  cent.      In  a   later  section,   this    is   shown  to   be 
important. 

Firing  of  the  high  voltage  pulser  causes    an  electrostatic    inter- 
action between  the  two  discharges.      This  may   result   from  capacitive 
coupling,   or  from  the  change   in  field  configurations   near  electrodes 
as    the  pulses    are   initiated.      In  either  case  the  mechanism   is   unclear, 
causing  an   increase   in  the   signal  detected   in  one   instance  and   a  decrease 
in  another    (figure   2-6). 

Figure  2-6  shows  the  form  of  the  acoustic  waves  detected  by  the 
diode  for  pulser  voltages  of  2  and  4  kv.  This  figure  shows  that  the 
higher  voltage  pulse  produces    a   larger  electrostatic    interaction,    a 
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Figure  2-6.   Traces  (a)  and  (d)  are  1001  Hz  square  waves  used 
as  base  lines  for  calculation  of  acoustic  wave  velocity. 
Nominal  time  base  is  0.5  msec/division.   Traces  (b)  and  (c) 
are  the  acoustic  waves  observed  by  the  diode  for  pulser  voltages 
of  2  and  4  kilovolts.   Density  increase  is  indicated  by  a  down- 
ward deflection  of  the  beam.   Deflections  prior  to  2.5  msec  are 
due  to  electrostatic  interaction.   The  velocity  of  the  acoustic 
wave  in  (b)  is  440  m/sec,  and  for  (c)  is  448  m/sec,  if  acoustic 
waves  are  assumed  to  arise  at  the  moment  of  maximum  light  in- 
tensity in  the  exciting  discharge.   Density  variation  in  (b) 
is  0.34  per  cent,  and  in  (c)  is  0.71  per  cent. 
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larger  amplitude  acoustic  wave,  and  an  acoustic  wave  which  arrives  earlier. 
Propagation  velocities  of  these  two  waves  are  440  m/sec  and  456  m/sec 
respectively,  if  the  waves  are  assumed  to  originate  at  the  time  of  break- 
down of  the  gas.   The  latter  velocity  corresponds  to  the  sonic  velocity 
of  small  amplitude  waves  for  a  neutral  gas  temperature  of  24°C.   This 
apparent  difference  of  propagation  velocities  cannot  be  resolved  without 
a  second  diode  in  the  tube  or  some  other  device  to  allow  time  of  flight 
to  be  measured  within  the  tube.   If  both  waves  are  assumed  to  propagate 
with  the  same  velocity,  75  microseconds  are  required  for  the  formation 
of  the  acoustic  wave  by  the  smaller  amplitude  pulse. 

Damping  of  acoustic  waves  by  collisions  and  heat  transfer  to  the 
walls  predominates  over  collision  and  heat  conductivity  losses  in  the 

volume  for  the  pressure  range  of  this  experiment.   25  1  The  damping 

-2    1 
coefficient  is  of  the  order  of  2  x  10   cm  . 
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CHAPTER  III 
ACOUSTIC  WAVE  INTERACTIONS 

A.   Translation  of  the  Gas 

Acoustic  waves  were  propagated  across  a  short  discharge  in  another 
tube  in  order  to  show  the  translation  of  the  gas  caused  by  these  waves. 
In  figure  3-1,  static  light  intensity  in  relative  units  is  plotted  as  a 
function  of  distance  from  the  centerline  of  the  discharge.   The  shape  of 
this  curve  is  approximately  Gaussian.   Arrows  indicate  photomultiplier 
locations  for  traces  (b)  and  (c)  of  figure  3-2. 

Trace  (a)  of  figure  3-2  shows  the  current  in  the  exciting  discharge 
as  a  function  of  time.   The  trace  begins  with  a  decrease  of  current, 
which  causes  a  rarefaction  wave  to  propagate  due  to  the  removal  of  a 
heat  source.   In  trace  (b),  the  photomultiplier  is  located  toward  the 
exciting  discharge  from  the  centerline  of  the  observing  discharge.   This 
causes  a  brighter  area  to  be  moved  in  front  of  the  slit  as  the  acoustic 
wave  arrives,  while  trace  (c)  shows  opposite  behavior  of  the  light  in- 
tensity as  a  darker  area  is  moved  in  front  of  the  photomultiplier. 
Effects  are  reversed  in  both  traces  as  current  in  the  exciting  discharge 
is  increased  and  a  compression  wave  is  propagated. 

When  the  acoustic  waves  were  propagated  along  the  axis  of  the  dis- 
charge tube,  a  similar  effect  was  noted  in  the  standing  striations  of 
cathode  regions.   However  the  light  intensity  was  nearly  square  wave 
modulated  when  E4  was  the  cathode  and  the  exciting  discharge  current  was 
modulated  by  a  square  wave.   The  behavior  with  E3  as  cathode  was  similar 
to  that  observed  for  propagation  of  the  acoustic  waves  transverse  to  the 
discharge  axis,  and  did  not  show  this  square  wave  modulation.   Figures 
3-3  and  3-4  show  the  waveforms  observed  when  each  of  these  electrodes 
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Lgure    3-1.       Light     intensity    as    a    function   of    position    trans- 
verse    to    the    positive    column   of    an    observing   discharge. 
Arrows    indicate     location   of    the    phot omul t ip 1 ier    for    traits 
(b)    and     (c  ;    of    figure    3   2.       For    this    experiment  ,    acoust  ie 
waves    were    propagated    across    the    column    (from    right     to 
left     in    this   d  lagram)  . 


27 


Figure    3-2.       (a)    Current    variation    in    the    ex<  ii  Lng    dis- 
charge  as   ,i   Function   oi    time;    50   ma /division.       (b) 
Light     intensity   variation  on   the   side    of   the 
olumn   of    the  observing   discharge    toward    the   exciting 
discharge.       (c)    Light     intensity   variation   on    the    side 
of    the   positive   column  of    the   observing   discharge    away 
from    the  exciting    discharge.       Current     in    the    obser- 
discharge   =    3.5    ma.      Time   base    =    2  msec/division. 


50  msec 


20  msec 


2  msec 


Figure  3-3.   Light  intensity  variations  on  the  cathode 
(exciting  discharge)  side  of  the  first  standing  stria- 
tion  in  the  positive  column  of  tube  II.   E3  is  cathode, 
0.7  torr,  16  ma.   Time  per  division  is  shown  to  the 
right  of  each  trace.   Rarefaction  wave  arrives  first 
in  each  case.   All  vertical  scales  are  identical. 
Tube  II  is  the  tube  used  for  acoustic  experiments. 
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Figure  3-4.   Light  intensity  variations  in  the  head  of  the 
positive  column  as  a  function  of  time.   Time  scale:  5 
msec/division.   (a)  Current  variation  in  exciting  dis- 
charge, 20  ma/division. (b)  Anode  side  of  first  standing 
striation.(c)  Cathode  side  of  first  standing  striation. 
(d)  Negative  glow.   E4  is  cathode,  0.7  torr,  16  ma.   Rare- 
faction wave  propagated  first. 
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is  operated  as  cathode.   Figure  3-3  also  shows  that  approximately  150 
msec  is  required  for  light  intensity  to  return  to  the  original  value. 
Amplitude  of  the  compression  phase  of  the  wave  is  decreased  as  frequency 
of  the  modulating  square  wave  is  increased. 

No  satisfactory  explanation  has  been  found  for  the  nearly  square 
wave  modulation  of  light  intensity  in  the  cathode  regions  when  E4  is 
used  as  cathode.   This  occurred  only  when  current  of  the  exciting  dis- 
charge was  square  wave  modulated. 
B.   Interaction  in  a  Homogeneous  Positive  Column 

Time-space  photographs  are  extensively  used  in  this  section  to  pre- 
sent the  several  effects  which  were  observed  when  acoustic  waves  were 
propagated  through  the  discharge.   For  these  photographs,  electrode  E3 
is  always  located  at  the  upper  limit  of  the  trace,  electrode  E4  at  the 
lower.1  •   Travel  of  the  photomult  iplier  is  40  cm  in  all  cases,  with  inter - 
electrode  separation  39  cm,  and  approximately  equal  overlap  at  each  end 
of  the  carriage  travel. 

The  translation  of  gas  discussed  in  the  preceding  section  gives 
rise  to  large  variations  of  light  output  for  regions  of  the  discharge 
having  standing  striations.   A  standing  striation  is  always  located  at 
the  head  of  the  positive  column,  and  other  standing  striations  which  may 
be  present  are  damped  in  the  direction  of  the  anode.   Damping  is  in- 
creased as  current  is  reduced,  leaving  the  column  nearly  homogeneous 
over  3/4  of  its  length  for  currents  below  about  2  ma. 

Interaction  in  the  positive  column  is  first  examined  by  using 
acoustic  waves  generated  from  square  wave  modulation  of  the  exciting 
discharge  current,  and  secondly  with  the  larger  amplitude  acoustic  waves 
from  the  high  voltage  pulser.   In  the  first  case,  the  general  effect  of 
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compression  and  rarefaction  waves  upon  the  column  is  established,  while 
with  the  larger  pulse  an  additional  effect  is  observed. 

Figure  3-5  shows  anode  to  cathode  propagation  of  acoustic  waves 
produced  by  square  wave  modulation  of  the  exciting  discharge.   A  rare- 
faction wave  is  propagated  first  in  all  experiments  with  the  square  wave 
modulation  technique,  and  it  is  seen  that  this  rarefaction  causes  an  in- 
crease in  light  emitted  by  the  discharge  (recall  that  increases  of  light 
intensity  are  recorded  as  dark  areas  in  the  photograph).   The  compression 
wave  causes  a  decrease  in  the  intensity  of  light  emitted  by  the  dis- 
charge.  Waves  of  stratification  also  originate  from  the  cathode  region. 
These  will  be  discussed  in  another  section. 

If  the  average  energy  of  the  electrons  remained  constant  during  the 
compression,  the  increased  collision  frequency  (at  the  higher  pressure) 
would  cause  an  increase  of  ionization  rate  and  an  increase  of  the  light 
emitted  by  the  discharge.   Since  the  light  intensity  decreases  in  the 
compression  phase,  the  electron  energy  must  decrease.   Increased  elastic 
collision  or  decreased  electric  field  may  explain  this  decrease  of  elect- 
ron temperature.   These  will  be  examined  in  a  later  section. 

Figure  3-6  shows  anode  to  cathode  propagation  of  an  acoustic  wave 
produced  by  the  pulse  discharge.   In  this  figure,  the  path  of  the  acoustic 
wave  is  shown  as  a  bright  line  preceded  by  a  thin  dark  line.   The  bright 
line  indicates  a  decrease  of  light  intensity  in  the  compression  phase, 
and  is  caused  by  a  decrease  in  the  electron  temperature.   The  dark  line 
indicates  an  increase  of  light  intensity  at  the  leading  edge  of  the  wave. 
Two  mechanisms  can  cause  this  increase: 

(a)  compression  with  constant  electron  temperature  (phase  lag 
of  temperature  change) 


32 


Path  of 
rarefaction 
wave 


Anode 


40  cm 


Slow  wave 

of 
strat  if  ica 
t  ion 


Electric 

Path  of 

Slow 

interaction 

compression 

wave 

with 

wave 

of 

exciting 

strat if i 

discharge  / 

s>  cation 

Cathode 


Figure  3-5.   Time  and  space  development  of  light  intensity  in  a  1.3 
torr,  2.6  ma  neon  discharge  resulting  from  anode  to  cathode  prop- 
agation of  rarefaction  and  compression  waves  produced  by  square 
wave  modulation  of  the  exciting  discharge.   Lighter  areas  in 
photo  correspond  to  darker  areas  in  the  discharge.   Time  scale 
is  2.0  msec/division. 


33 


Path  of 
acoustic 
disturbance 


Fast    wave 
of 
stratification 


Anode 


Cathode 


Slow  waves 
of 
stratification 


Figure   3-6.      Time-space  development   of   light   intensity   in  a  2.2 
torr  neon  discharge  at   3.1  ma.      Anode  to  cathode  propagation 
of   a   compressive   acoustic  disturbance.      Time  scale   is   2  msec/ 
division.      Thin  dark   line  preceding  bright   line   in  acoustic 
interaction   indicates   a  slight    increase  of   ionization  rate. 
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Figure   3-7.     Time-space  development  of   light    intensity   in  a  0.82 

torr   35  ma  neon  discharge.      Anode  to  cathode   sound   propagation. 
Time  scale   1  msec/division.       Striation  speeds   are   140m/sec, 
66  m/sec,    33  m/sec   respectively.      The  fastest    striations   appear 
to   be  excited   by  travel   of   the  acoustic  wave   along  the  positive 
column. 
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(b)  increase  of  electric  field  and  electron  temperature. 
This  phenomenon  was  not  observed  in  all  photographs  with  the  high-voltage 
pulser,  and  was  never  seen  when  square  wave  modulation  of  the  discharge 
current  was  used  for  acoustic  wave  production. 

Anode  to  cathode  propagation  of  the  acoustic  wave  causes  changes  of 
ionization  rate  in  time  and  space  which  are  similar  to  those  present  in 
moving  striations.   Figure  3-7  shows  striations  which  are  apparently 
produced  during  propagation  of  the  acoustic  waves  along  the  column.   By- 
viewing  this  photograph  along  the  apparent  path  of  the  acoustic  wave, 
with  the  eye  near  the  plane  of  the  page,  a  series  of  lines  can  be  seen. 
The  acoustic  wave  arrives  at  the  anode  1.6  msec  after  the  beginning  of 
the  sweep.   At  this  time,  light  intensity  is  observed  to  increase  (dark 
area  in  photograph).   This  increase  can  be  followed,  with  some  inter- 
ruptions, for  most  of  the  distance  of  the  column.   A  second  dark  area 
parallel  to  this  is  also  visible  at  the  anode  end  of  the  photograph. 
The  position  at  which  the  third  dark  line  should  appear  is  somewhat  ob- 
scured by  the  bright  areas  associated  with  the  maximum  density  portion 
of  the  acoustic  wave.   Succeeding  striation  patterns  become  clearer. 

Striation  production  by  the  acoustic  waves  in  an  inhomogeneous 
positive  column  is  considered  in  the  following  section. 
C.   Interaction  in  an  Inhomogeneous  Positive  Column 

As  the  acoustic  wave  encounters  standing  striations  in  the  positive 
column,  changes  of  the  ionization  rate  at  these  points  also  produce 
moving  striations.   Light  intensity  gradients  are  small  for  all  but  the 
first  few  standing  striations,  so  that  translation  of  the  gas  will  pro- 
duce negligible  variation  of  the  photomult ipler  signal  in  that  half  of 
the  positive  column  toward  the  anode.   Figure  3-8  shows  the  standing 
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Figure  3-8.   Light  intensity  as  a  function  of  position  in  a  0.82 
torr,  15  ma  neon  discharge.   Standing  striation  patterns  had 
similar  form  at  10  ma,  smaller  amplitudes.   Arrows  show  approx- 
imate phases  of  standing  striations  at  which  moving  striations 
of  figure  3-9  originate.   (Small  variation  as  light  decreases 
to  zero  at  anode  end  is  due  to  a  layer  of  sputtered  material 
on  the  glass. ) 
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striation  configuration  in  a  15  ma  discharge  at  a  pressure  of  0.82  torr. 
Figure  3-9  shows  the  production  of  moving  striations  by  the  acoustic 
wave  as  it  encounters  the  standing  striations.   These  decay  and  merge 
directly  into  the  wave  of  stratification  which  is  excited  from  the  stand- 
ing striations  nearer  the  cathode.   Wavelength  of  the  standing  striations 
is  twice  that  of  the  moving  striations. 

Amplification  of  the  wave  of  stratification  increases  as  discharge 
current  is  brought  near  a  critical  current  (at  which  self -excited  stria- 
tions occur).   Waves  of  stratification  can  be  excited  by  smaller  local- 
ized changes  of  ionization  rate  in  these  cases.   Figure  3-6  shows  pro- 
duction of  a  fast  and  a  slow  wave  of  stratification  from  the  region  near 
the  second  standing  striation,  and  a  second  slow  wave  from  the  cathode 
region.   In  both  cases  it  is  difficult  to  fix  the  point  of  origin  of  the 
wave.   A  third  slow  wave  is  visible  at  the  right  hand  side  of  the  photo- 
graph.  This  is  due  to  oscillation  of  tube  voltage  and  current  as  the 
first  two  slow  waves  arrive  in  the  anode  region.   Figure  3-7  also  shows 
waves  of  stratification  originating  from  more  than  one  region  of  the 
discharge. 

At  lower  currents,  amplitude  of  the  standing  striations  is  smaller, 
and  waves  of  stratification  originate  only  in  the  cathode  region.   Figure 
3-5  shows  that  waves  of  stratification  originating  in  the  cathode  region 
due  to  compression  and  rarefaction  acoustic  waves  have  nearly  opposite 
courses  of  development  of  light  intensity.   The  change  of  electric  field 
in  the  cathode  region  causes  the  wave  of  stratification  to  be  excited. 
Since  the  compression  and  rarefaction  waves  have  opposite  effects  on  both 
of  these  quantities,  the  waves  of  stratification  which  are  produced  should 
be  opposite  in  their  course.   Theory  and  experiment  have  shown  that  the 
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Path  of  acoustic 
wave 


Slow  wave  of 
stratification 


Anode 


Cathode 


Figure   3-9.      Time-space  development  of   light   intensity   in  a  0.82 
torr,    10  ma  discharge.      Anode  to  cathode  propagation  of  the 
acoustic  wave.      Time  scale  is   1  msec/division.      Striation 
production  by  sound  wave   is  evident. 
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Fast   and   slow  waves   of 
stratification   from 
electric   action 


Anode 


Slow  waves  of 
stratification  from 
acoustic  interaction 


Cathode 


Path  of 


rarefaction  wave 


Path  of 
compression  wave 


Figure  3-10.  Time-space  development  of  light  intensity  in  a  0.7 
torr  16  ma  neon  discharge.  Anode  to  cathode  propagation  of 
acoustic  disturbances.   Time  scale  is  2  msec/division. 


HO 


wave  of  stratification  is  an  "odd"  (signed)  effect. |  43,  48 

The  effect  of  direction  of  propagation  is  illustrated  by  comparing 
figures  3-10  and  3-11.   In  figure  3-10,  the  wave  of  stratification  pro- 
duced in  the  cathode  region  following  anode  to  cathode  propagation  of 
the  acoustic  wave  is  highly  damped.   For  cathode  to  anode  propagation  of 
the  acoustic  wave,  striations  appear  along  its  path.   This  behavior  was 
also  noted  in  a  4.0  torr  discharge,  where  striations  associated  with  the 
fast  wave  of  stratification  were  excited  along  the  path  of  the  acoustic 
wave  (figure  3-12).   Group  velocity  of  the  wave  of  stratification  ia  ; 
160  m/sec,  significantly  smaller  than  the  sonic  velocity,  and  the  stria- 
tions excited  by  the  acoustic  wave  are  attenuated  before  being  amplified 
again  in  the  wave  of  stratification. 

In  figure  3-13,  the  acoustic  wave  is  shown  to  decrease  the  ampli- 
tude of  self-excited  fast  striations.   These  striations  are  associated 
with  direct  ionization  of  the  neutral  gas  molecules,  and  the  decrease  of 
electron  temperature  caused  by  the  acoustic  wave  reduces  the  amplitude  of 
the  striations. 
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Path  of 
rarefaction  wave 


Path  of 
compression  wave 


Cathode 


Anode 


Fast  and  slow  waves  of 
stratification  from 
electric  interaction 


Slow  waves  of 
stratification  from 
acoustic  interaction 


Figure  3-11.   Time-space  development  of  light  intensity  in  a  0.7 
torr,  16  ma  neon  discharge.   Cathode  to  anode  propagation  of 
acoustic  disturbance.   Time  scale  is  2  msec/division. 
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Cathode 


40  cm 


Path  of  acoustic  wave 

\ 


Anode 


Fast  wave  of    — 
stratification 


Slow  wave  of 
stratification 


Figure   3-12.      Time-space  development   of   light   intensity   in  a  4.0 
torr,   0.2  ma  neon  discharge.      Cathode  to   anode  sound   propa- 
gation.     Time   scale:    1  msec/division. 
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Time  of  arrival  of 
acoustic  wave 


Figure  3-13.   Light  intensity  variation  with  time  at  a  distance  of 
1.3  cm  in  front  of  the  anode.   Discharge  currents  (a)  2.2  ma 
(b)  2.1  ma;  pressure  0.77  torr.   Cathode  to  anode  sound  propa- 
gation.  Time  scale  2  msec/division,  light  intensity  traces  to 
the  same  relative  scale.   Light  intensity  decreases  as  compressive 
acoustic  wave  arrives. 
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CHAPTER  IV 
ANALYSIS  OF  ACOUSTIC  INTERACTIONS  IN  THE  DISCHARGE 

A.   Change  of  Electron  Temperature 

The  experiments  showed  that  the  light  emitted  by  the  homogeneous 
positive  column  decreases  during  the  compression  phase  of  the  acoustic 
wave,  and  that  this  must  be  explained  on  the  basis  of  a  change  in  temper- 
ature.  In  some  cases  the  light  emitted  by  the  column  increases  for  a 
short  time  preceding  this  decrease. 

Since  the  ions  are  nearly  in  thermal  equilibrium  with  the  neutral 
gas,  neutral  gas  density  perturbations  will  tend  to  produce  the  same 
fractional  changes  in  ion  density.   These  can  be  expressed  as 

N._ 
n+  =  —  •  nn  (11) 

Nno 
where 

nj  =  nj-njo  h|"Nj0  •  (i2> 

3 
Frequencies  of  the   acoustic  waves   were  <—  10     Hz,    giving   a  charaoter- 

o 

istic  time  of  <~  10        sec   for  the   ion  density  perturbation.      Ambipolar 
diffusion  times   are  of  the  order  of 

T~  ^-£ —     =  4  x  10-5   sec    ,  (13) 

a  '*■    D 
a 

indicating  that  ion  density  changes  will  be  removed  through  ambipolar 
diffusion  to  the  walls.   However,  Saxton  observed  ion  density  variations 
of  the  form  given  by  equation  (11).   Diffusion  of  the  charged  particle 
density  perturbation  evidently  does  not  take  place. 

The  increase  of  neutral  particle  density  also  decreases  the  mobil- 
ity of  the  charged  particles,  and  increases  electron  energy  losses  due 
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to  an  increase  of  the  collision  rate.   These  tend  to  bring  about  a  de- 
crease in  electron  temperature.   Electric  field  may  also  be  varied  by 
the  density  change,  and  thus  lead  to  an  additional  change  in  electron 
temperature.   In  most  cases,  electron  temperature  is  related  to  the  re- 
duced electric  field  (E/p) ,  so  that  it  is  difficult  to  separate  the 
effect  of  field  and  pressure  variations  on  the  electron  temperature  for 
this  experiment. 

Using  the  relation  between  electron  temperature,  pressure  (density) 
and  electric  field  given  in   17  I,  it  is  possible  to  examine  the  effects 
of  neutral  gas  density  variations.   This  equation  is 

Te  =  -*4-  -%=     .  E  (14) 

where 

Te   is  in  eV 

^   electronic  mean  free  path 
=  1/<T-  Nn 

f_   fraction  of  electron  energy  lost  per  collision. 
This  equation  can  be  written 

T„  =  -JET.  ■   1 t-  (15) 


3  •  <?"   /t> 


Von  Engel  graphs  the  variation  of  f_  as  a  function  of  E/p  in   17  . 

The  curve  is  for  a  constant  neutral  gas  temperature  and  is  therefore 

equivalent  to  f_  as  a  function  of  E/Nn.   For  pressures  and  fields  in  the 
experimental  range  this  dependence  can  be  expressed  as 

f   =  f_0  (— )b  .  (16) 

P 

Variation  of  electron  temperature  with  neutral  gas  density  can  therefore 
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be   determined   from 


In  general 


and 


so 


dT 


e 


6Te 


<*Nn  d(E/p) 


6% 


6  (E/p) 


d(E/p) 


dCETp) 


dNn 


<*Te      _      Te 


3/Ttf"   /iT 


j_  r  ^e 


~5X 


[-+H 


— ] 

Nn  J 


(17) 


(18) 


(19) 


(20) 


In  the   succeeding   paragraphs,    it   will   be   shown  that     <jE/dNn-0.      With 
this   knowledge,    the  temperature   change  with  density  depends   upon  b  as 
follows 


b  <   2 


dTc 


6k 


*~  o 


b   =   2 


b  >  2 


dT. 


dNf 


=   0 


-7  0    . 


For  the  case  of  neon  in  the  experimental  range  of  reduced  field  (0.8  - 
5.0  v/cm  •  torr) ,  b  =  1.9.  Thus  the  variation  of  electron  temperature 
with  pressure   is   small. 

Electric   field  variation  can  be   analyzed  on  the  basis   of  the  work 
of  Pekarek   and  Masek.    1 5 8  I      In  their  method,    the   electron  density  varia- 
tion  is   assumed  to   follow  the   ion  density  variation  quasi-statically, 
with   its   form  governed   by  the  continuity  equation: 
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dN   +  v  (NV)=  Q.  (21) 

When  perturbations  of  the  form 

n_   =   N_  -  N_Q      |n_|  ^   N_Q  (22) 


e   =   E  -  EQ       |e  |  «  EQ  (23) 


are  substituted,   this   equation  becomes 


d2n_  dn_  ^ 

_ +    j  .     T7         -4-     /y        NT  O 


e 


0    =   D_  —  +>w   En        r-      +//     N0   -^ 

£z2       ^"  °       <^Z        ^-      °  ^Z 


+  N    E      fc*-      ,-^IL-  . 

°   °  ^       <*Z 


(24) 


Electric   field   is   related   to  charged   particle  density  through  Poisson's 
equation 

4s-  =     —   K   -  n->    •  C25) 

0Z  4?T  CQ 


Since 


M     =  -£—      ,  (26) 

N„ 


then 


and  equation    (24)    can  be  written    in  the  form 


(27) 


Using  the  two-sided  Laplace  transform  defined  in  61 I  by 


(28) 


-  oo 
the  equation  for  the  electric  field  is 


f(s)  =  s  /   h(t)  e~stdt  ,  (29) 
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e   = 


^TT€„ 


1    - 


1    -    ^s2    J 


(30) 


where   e   is   the  operator    image  of  e.      Equation   (30)    is   to  be  compared 
with  Pekarek  and  Masek's   equation    (10)    which   is 


e   = 


J_  (_2o__)  h L_± 


where    j^   =   2D_///_EQ, 


(31) 


(32) 


The   acoustic   perturbation  therefore   gives   a  result   similar  to    letting 
f%    -*■  oo,    for  which  the  change  of  electric   field   is   zero. 

If  the  change  of  mobility  is   not  taken   into   account,    the  compressive 
acoustic   perturbation  can  be  handled  through  their  calculation  for 
/,    =   100  ^D.      This  causes   an  electric   field  change  of   -2.5   x  10~3  v/cm 
for  a   fractional  density  variation  of  0.5   per  cent. 

Comparing  the  terms    in  the    last   bracket   of   equation    (20),    and 
assuming  E      =   2  v/cm 

Nn  =    3.5  x  1016   cm"3, 
we   find 

=    -    1. 3  x    10 


AN, 


J-  =5.7   x    10"17 


Thus  the  terms   are   of  the  same  order  of  magnitude.      It   has   been  assumed 
that  the  perturbation  of    ion  density  does   not  decay  due  to   ambipolar 
diffusion  to   the  walls.      If   ambipolar   diffusion  does   prevent  this   field 
change,   there    is   still   a  temperature  variation  resulting  from  the  neu- 
tral  gas   density   perturbation.      Such   a  change  must    be   attributed  to 
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variation  of  collision  losses  as  the  collision  frequency  is  changed. 

Magnitude  of  the  temperature  change  can  be  obtained  from  equation 
(20).   Recasting  it  in  the  form 


ATe  b    I  AE     AN, 


■  ■  4-H- 


Te  E      Nn 


(33) 


of 
the  effects  fractional  changes  of  electric  field  and  neutral  particle 

density  are  readily  seen.   For  the  conditions  given  in  the  preceding 

paragraph 

t  Te   =  -9.3  x  10"4  .  (34) 

Te 

The  predicted  temperature  variation  is  therefore  of  the  order  of  one 

3 
part  in  10  .   Sensitivity  of  the  time  space  display  method  is  quite  evi- 
dent from  this  value. 

The  short  duration  increase  of  light  intensity  which  sometimes  pre- 
ceded the  longer  decrease  in  the  compression  phase  may  be  due  to  the 
effect  of  the  finite  length  of  the  acoustic  waves.   This  finite  length 
wave  can  be  described  by 

f  =  -f  _  sin(kZ  -cut)   U(t —  )  (35) 

for  a  wave   which   is   propagating   from  anode  to   cathode  with  velocity  c, 
and   a  compression  phase   at    the  leading  edge.      In  this  equation 

J    =  Nn  "Nno    ■  |Sh<Nno  (36) 

and  the  wave   is  of   semi-infinite   extent.      An   exponential  damping  term 
should   be   included   to   give  a  form  similar  to  that   shown   in  figure   2-5. 
The   change  of   ionization  rate  can  be   expressed   as 


dz      _  dz  dNn 


dt  dNn  dt 
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(37) 


and  from 


M 


dz         3/2  _ 


=  -  c,x   e~*  (38) 


dNn 


where 


Thus 


x  =  eV,-/kT 

Am  C 


4^—   =  c,  %     x3/2e~X   -cocos(kZ  -  cot)  U(t  -  -) 
at     i  ->o  v 


[ 


-4 


+  sin  (kZ  -  cut)  •  S(t  -  2)  ,  (39) 


and  the  solution  will  have  a  transient  which  propagates.   This  transient 
may  explain  the  slight  increase  of  ionization  rate  sometimes  observed 
at  the  leading  edge  of  the  acoustic  wave. 
B.   Excitation  of  Acoustic  Waves  by  the  Discharge 

Ingard  has  predicted  that  acoustic  waves  can  be  generated  in  the 
neutral  gas  of  a  discharge  through  collisions  with  the  charged  particles, 

27   Working  from  the  linearized  equations  for  continuity  of  mass, 
momentum  and  energy  of  the  neutral  gas,  with  charged  particle  inter- 
actions entering  as  source  terms,  he  obtains  the  equation 

-4-  -S?E_  -  V2  P  =  £$-  4s-  -  V-  F  .  (40) 

^2-  6t2  c2  ^r 


In  this   equation 


c        sonic  velocity 

=  VYkTn/Mn  (41) 

p  variation   in  pressure 

/  ratio  of   specific   heats 

H  heat  energy  source   term 

momentum  transfer   source  term. 
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This  equation  shows  that  acoustic  waves  can  be  excited  by  a  change  of 
the  heat  transfer  rate  or  a  spatial  gradient  in  momemtum  transfer. 
Ingard  attributed  acoustic  disturbances  from  moving  striations  to  the 
former,  while  Briand  and  Delcroix  observed  acoustic  waves  which  they 
ascribed  to  momentum  transfer  from  the  ions  to  the  neutral  gas.  [  7 

Before  carrying  out  a  calculation  to  compare  the  magnitudes  of 
these  two  effects  in  the  positive  column  when  moving  striations  are  pre- 
sent, it  should  be  pointed  out  that  variation  of  the  power  dissipated  in 
the  cathode  fall  due  to  a  change  of  current  will  cause  a  much  larger 
acoustic  wave  than  any  change  in  the  conditions  of  the  positive  column. 
Normal  cathode  fall  for  neon  with  a  molybdenum  electrode  is  approximately 
110  volts;  I  17   the  cathode  properties  of  molybdenum  and  tantalum  are 
similar.   Potential  difference  between  the  first  two  standing  striations 
in  the  neon  positive  column  is  18.5  volts.  I  16   This  is  larger  than  the 
potential  difference  between  peaks  of  moving  striations.   Assuming  all 
energy  dissipated  in  each  region  is  consumed  in  gas  heating,  a  step- 
wise change  in  discharge  current  will  cause  an  acoustic  wave  to  origi- 
nate from  the  cathode  region  which  will  be  larger  by  a  factor  of  5  than 
that  which  will  be  produced  from  the  region  between  the  standing  stria- 
tions.  However,  it  is  possible  to  have  moving  striations  present  in  the 
discharge  with  no  variation  in  discharge  current,  or  to  have  current 
variation  small  enough  that  the  acoustic  waves  produced  in  the  cathode 
region  are  smaller  than  those  which  arise  from  the  moving  striations. 


[lO,  ll] 


Consider  the  case  of  acoustic  waves  produced  by  perturbations  of 
an  infinite,  homogeneous  dc!  discharge  column.   In  order  to  compare  the 
magnitudes  of  the  source  terms  in  the  equation  derived  by  Ingard,  it  is 
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assumed  that  the  electric  field  is  2.0  volts/cm,  electron  temperature 

is  2.0  eV,  and  neutral  gas  pressure  is  1.1  torr.   Under  these  conditions, 

the  following  parameters  are  applicable  J  17 

N_  =  2.4  x  1016  m-3 
aa^   =  1.2  m2/volt  sec 
yU.   =  80  m2/volt  sec 


\   =  1.2  x  10 


-4 


m 


\   =    6.6   x   10~4  m   . 
As   an   approximation  to  the  momentum  transfer  term,    it   can  be  as- 
sumed that   all  momentum  gained  by   a   charged   particle  from  the  electric 
field   in  traveling  one  mean  free  path  along  the  direction  of   the 
electric   field   is  given  to   a  neutral   particle  at   the  end  of  that   free 
path.      This    corresponds   to   the   fact   that  distribution  functions  of 
charged   particles    in  electric   fields   can  be  divided   into   a  part  which   is 
symmetric    (isotropic)   and   a  part   which   is   asymmetric.      The   average  energy 
of   the   symmetric   part    is   much   larger,    and  for  electrons,   corresponds  to 
the  temperature  measurable   by  an  electric   probe.      Average  velocity  of 
the  asymmetric   part   corresponds  to   the  mobility  or  drift  velocity.      As- 
suming that   all   particles   have  the  drift  velocity  of   the  particular 
species,    the  momentum  transfer   is   given  by 

"F   =   n(m+    .   v+    •   v>c+  +  m_   •  "v_    •   yc   ).  (43) 

The  effective  collision  frequency   for  asymmetric  momentum  transfer   is 


>b+   =        |vj     /Xf      >  (44) 


and   as   the  elastic   collision  cross-section  of  neon   is   nearly   independent 
of   energy,   the  mean   free  path   is   essentially  constant.      The  velocities 
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used   in  calculating  these  collision   frequencies    are  mobility  velocities 

v+   =      ±//+  E  (45) 

since  only  collisions  along  the  direction  of  drift  motion  are  to  be 
considered  for  momentum  transfer.   This  gives 

F  =  n(  \*~   -   m-^  )  E  .  |  E 1  (46) 

Comparison  of  the  magnitudes   of   ion  and  electron  contributions   to  this 
equation  shows   that 

m.    *  M.    '    3 

^    45,  (47) 


m 


T 


and  the  electrons  can  be  neglected.   Note  that  this  is  independent  of 
electron  and  ion  temperature.   Thus  the  equation  of  interest  is 

V-    F  =   2  .  n  .  f .   E   .  -AlL-  .  (48) 

A+  I   I     AX 

Large  changes   of   electric   field   are  observed    in  the   presence  of  moving 
striations,    and  an  assumption  of 


— —     =      1   volt/cm 
Ax 

is  well  within  the  range  found  in  other  experiments.   12,  67    With 


these  values 


V-  F  =   38.4  joules/m5  (49) 


To  compare  this  with  the  heat  energy  term  resulting  from  the 
change  of  electron  temperature  in  moving  striations,  it  is  assumed  that 
electron  temperature  is  doubled  in  20  microseconds,  that  each  electron 
has  the  energy  of  the  distribution  function  average  and  that  all  electrons 
collide  elastically  with  neutral  atoms.   The  fraction  of  energy  given  up 


54 


by  each  species  in  an  elastic  collision  is  denoted  f+.   Relative  con- 
tributions of  ions  and  electrons  can  be  approximated  by 

H_  n_    •    f_    •      Vc_    .    kT_  ,cn. 

=  26.  (50) 


H+  n+    •    f+    •      yC+    •    kT+ 

Thus   the   ion  contribution  can  be   neglected,    especially  since  the  tempera- 
ture of  the  electrons  doubles  while  that  of   the    ions  changes    little. 

Ingard   gives   the  following   equation  to   express  the  heat   contribu- 
tion of   the  electrons    [ 2 7 


/-i     dn    -  P 

~   ro 


c 


2       dt 


m0 
(-^-) 


^lJfS3\^l.  (51) 


Assuming  the  temperature  change  to  be  the  dominant  factor,  this  equa- 
tion becomes 


c 


where 


<<*"'>  =    6    •     (  /  -    DC-l-y^)    =   5.2     <<T> 
«3->  =    3.95   x    1(T20  m2 
Ne      =    2.4   x    1016   m~3 
c        =456  m/sec 
Tn      =    300   °K. 
Using  the  values,    and  the  postulated  temperature  change, 

7-1      _^H_     _      3>52   x    1Q2    joules/m5.  (53) 

c2       C)t 

This    is   an  order  of  magnitude   larger  than  the   contribution  from  the 

momentum  term.      Discharge  conditions   under  which  E,     **"...     and   sftt 

°  Ax 

are  larger  and  A Te  is  smaller,  will  change  these  comparative  magnitudes, 
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The  experiment  of  Briand  and  Delcroix  was  performed  with  a  plasma 
produced  by  an  rf .'  discharge.   The  voltage  between  two  grids  bounding 
this  plasma  was  varied  at  a  given  frequency,  and  a  microphone  located 
normal  to  the  plane  of  these  grids.   Acoustic  waves  were  detected  at 
the  frequency  of  this  modulating  voltage.   Since  the  rate  of  heat  gen- 
eration is  unchanged  by  the  small  modulation  signal,  these  acoustic  waves 
must  be  generated  by  the  term  V*  F,  i.e.,  a  net  transfer  of  momentum  to 
the  neutral  gas  by  the  ions.   Amplitude  of  these  waves  was  observed  to 
increase  as  the  R-F  pow6r  increased.   Charged  particle  density  increases 
as  this  power  increases,  and  as  V-    F  is  linearly  dependent  upon  this 
density,  the  experiment  and  theory  are  in  agreement. 
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APPENDIX  I 
DISCHARGE  PROCESSES 

This  appendix  is  included  for  the  reader  who  is  unfamiliar  with 
glow  discharges.   The  discharge  is  described  generally,  and  then  dis- 
cussed from  a  microphysical  point  of  view.   Particular  attention  is  paid 
to  processes  important  to  wave  of  stratification  theory.   Cobine  |13j, 
von  Engel   17  and  Loeb  I 37 J  all  treat  the  discharge  processes  in  more 
detail,  and  the  stratification  work  has  recently  been  reviewed  by  Oleson 
and  Cooper  I  40 I . 
A.   General  Description 

Figure  4-1  shows  schematically  the  major  areas  of  a  low-pressure, 
cold  cathode  dc  glow  discharge.   A  discharge  of  this  type  in  any  vessel 
has  a  cathode  fall  of  definite  length  determined  by  the  gas  pressure  and 
type,  and  by  the  electrode  material.   This  region  is  essential  to  the 
production  of  electrons  for  the  discharge.   The  positive  column  is  the 
next  important  region.   A  glow  discharge  can  be  operated  without  a  posi- 
tive column.   Whenever  the  discharge  is  longer  than  the  cathode  fall, 
the  positive  column  fills  the  remaining  space  between  the  negative  glow 
and  the  anode.   It  is  this  part  of  the  discharge  which  is  used  for  plasma 
research  when  a  weakly  ionized  medium  is  suitable.   Commercial  applica- 
tions of  the  positive  column  use  its  light  output  in  "neon"  signs  and 
fluorescent  lights.   The  anode  fall  region  produces  ions  for  the  positive 
column,  and  is  not  important  unless  this  fall  of  potential  becomes  un- 
stable.  In  this  case,  anode  spot  oscillations  develop  and  disturb  the 
properties  of  the  positive  column. 

Current  in  the  discharge  can  be  varied  over  more  than  three  decades 
(0.5  to  500  mA)  with  constant  voltage  drop  across  the  discharge  tube. 
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Figure  4-1.   Characteristic  light  and  dark  areas  of  a  low  pressure 
dc  glow  discharge.   Lower  sketch  shows  electric  field  as  a 
function  of  position  in  the  positive  column,  with  the  coordi- 
nate system  chosen  such  that  the  positive  X  direction  is  in 
the  direction  of  the  electric  field. 
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This  current-voltage  behavior  characterizes  the  normal  glow  regime,  and 
finds  application  in  voltage  regulator  tubes.   For  discharge  currents  of 
several  hundred  milliamperes ,  heating  of  a  cylindrical  cathode  becomes 
noticeable.   Heating  of  the  tube  walls  also  becomes  appreciable  at  these 
higher  currents,  particularly  in  narrow  tubes.   Sputtering  of  cathode 
material  onto  the  walls  of  the  discharge  tube  presents  a  problem,  however, 
particularly  as  pressure  is  reduced  below  two  torr. 

The  light  emitted  by  each  section  of  the  discharge  is  characteristic 
of  that  region,  and  increases  with  increasing  discharge  current.   Light 
from  the  negative  glow  in  neon  is  always  orange,  as  is  that  from  the 
anode  fall  when  this  is  visible.   In  the  positive  column,  a  definite 
change  of  color  from  orange  to  red  is  observed  as  either  current  or 
pressure  is  increased.   A  standing  striation  always  heads  the  positive 
column,  and  several  other  spherically-shaped  light  regions  can  be  ob- 
served in  the  column  at  lower  pressures  within  the  experimental  range. 
These  are  highly  damped  toward  the  anode.   Viewing  the  striations  with 
a  prism  spectrograph  shows  that  lower  energy  (red)  lines  are  excited 
closest  to  the  cathode  with  orange  and  green  lines  appearing  toward  the 
anode. 

Thermionic  emission  of  the  electrons  from  a  heated  tungsten  wire  or 
heated  oxide-coated  cathode  eliminates  the  cathode  fall  region.   Other 
surfaces  will  produce  electrons  without  heating  and  allow  a  low  potential 
difference  to  be  used  in  maintaining  the  discharge.   24   Cold  cathode 
operation  was  chosen  for  these  experiments  for  several  reasons:  (1)  no 
external  power  supply  is  needed  as  in  the  case  of  the  heated  cathodes; 
(2)  electrodes  may  be  identically  constructed  and  the  polarity  reversed 
to  check  for  any  directional  dependence  in  the  tube;  (3)  simplicity  of 
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construction;  (4)  insensitivity  of  this  cathode  to  any  environment  as 
opposed  to  the  oxide  coated  cathodes. 

Regardless  of  choice  of  cathode,  the  negative  glow  is  always  pre- 
sent and  it  is  this  glow  that  gives  the  discharge  its  name.   Approximately 
equal  numbers  of  positive  and  negative  charges  are  present,  making  this 
the  first  region  with  plasma  characteristics. 

The  plasma  of  the  positive  column  is  a  slightly  ionized  gas,  with 

-6      -4 
fraction  of  ionization  10   to  10   .   Collisions  with  neutral  atoms 

play  a  dominant  role  in  this  region.   Mobility  and  ambipolar  diffusion 
(described  later)  govern  the  motion  of  charged  particles  in  the  column, 
and  recombination  is  assumed  to  be  significant  only  at  the  walls.   A 
small  positive  gradient  of  potential  in  the  direction  of  anode  is  essen- 
tial for  the  production  of  the  charged  particles  which  maintain  the 
discharge,  and  it  is  this  positive  gradient  which  names  this  portion  of 
the  discharge. 
B.   Cathode  Fall 

Electron  production  at  the  cathode  is  accomplished  primarily  through 
acceleration  of  ions  into  the  metal  surface  by  the  high  electric  field 
present  in  this  region  (figure  4-1).   As  the  ejected  electrons  are  accel- 
erated away  from  the  surface  by  this  same  field,  they  excite  neutral 
atoms  and  form  the  cathode  layer.   Further  acceleration  raises  electron 
energy  above  the  ionization  Level  and  leads  to  production  of  more 
electrons,  this  time  in  the  Crookes  dark  space.   These  electrons  repeat 
this  chain  of  events,  and  the  multiplication  process  grows  exponentially 
with  distance  from  the  cathode. 

If  we  let  ex  denote  the  number  of  ionizations  per  unit  length  of 
travel  of  an  electron,  and  assume  that  the  ionization  per  length  is 
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proportional  to  concentration,  then  a  single  electron  emitted  at  the 
cathode  will  give 


d 
(f o<  dx)  (54) 

o 


a 

f 

•*  n 


electrons   at    a  distance  d  from  the  cathode.      Assuming  also,    that    !/•/ 
ions   must    impinge  upon  the  cathode   to  cause  emission  of  a  single  electron, 
in  the   steady   state  each  electron  emitted  must   create   !/•/     ion  pairs   to 
maintain  the  discharge.      Thus 

d 
yfexp(j   c*  dx)    -   l]      =   1  (55) 

and 

<*  dx  =   ln(l    +      1   )  .  (56) 
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The   integral   gives   the   number  of   ionizing  collisions   from  the  cathode 

to   d.      This   can  also   be  approximated   by 

d 
<x    dx   =    >|Vn    ,  (57) 


a 


where    h    is   the  number  of    ion  pairs    per  electron  per  volt   of   potential, 
and   V     is   the  normal   cathode  fall: 

Vn  =  4-     ln(1  +  -*->• 

f\  1  (58) 

This  equation  shows  that  normal  cathode  fall  increases  as  ionization 
efficiency  decreases,  and  as  the  efficiency  of  production  of  secondary 
electrons  decreases.   For  the  case  of  neon  gas  and  an  aluminum  cathode, 
the  normal  cathode  fall  is  120  volts,  and  the  length  of  the  cathode  fall 
at  one  torr  is  0.64  cm.   16   Tantalum  has  a  slightly  lower  fall  of  poten- 
tial and  shorter  dark  space,  as  its  characteristics  are  similar  to  that 
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of  molybdenum.   3   Complete  data  on  tantalum  were  not  available. 

Cathode  potential  fall  can  be  reduced  and  current  density  increased 
by  improving  the  secondary  emission  coefficient — the  term  /  in  equation 
(58).   Although  ions  have  been  mentioned  as  producing  the  electrons, 
ultraviolet  photons  can  also  produce  electrons  through  the  photoelectric 
effect.   This  effect  is  enhanced  by  using  two  large  closely-spaced 
parallel  plates  as  cathodes,  and  is  called  the  hollow  cathode  effect. 
It  may  also  be  obtained  by  using  a  cylindrical  shell  of  small  diameter, 
and  it  was  in  this  manner  that  the  hollow  cathode  effect  was  used  in 
these  experiments.   Investigations  which  have  been  reported  show  the 
importance  of  separation  distance  for  flat  plates,  and  the  variation  of 
this  distance  with  the  gas  and  the  pressure.   78,  79   Although  no  data 
was  available  for  the  cylindrical  shells,  the  flat  plate  data  gave  an 
idea  of  the  minimum  and  maximum  diameters  which  could  be  used  in  the 
pressure  range  of  these  experiments. 

Magnetic  fields  placed  such  that  the  magnetic  lines  are  parallel 
to  the  plane  of  the  electrodes  have  been  used  in  conjunction  with  the 
hollow  cathode  construction.  |62,  63   These  fields  have  the  same  effect 
as  increasing  pressure,  provided  certain  limits  are  observed.   Since  an 
equivalent  pressure  increase  can  be  obtained,  the  effective  lower  pressure 
limit  for  a  given  geometry  can  be  reduced. 

As  electrons  leave  the  high  field  region  of  the  cathode  fall,  their 
average  energy  is  reduced.  The  light  emitted  in  the  negative  glow  shows 
this,  as  higher  energy  spectral  lines  appear  nearest  the  cathode. 

Sputtering  of  cathode  material  onto  the  discharge  tube  walls  hinders 
observation  of  cathode  area  light  intensity.   It  has  a  much  more  serious 
effect  in  the  case  of  a  sealed  tube,  for  the  pressure  is  decreased  as  gas 
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molecules    are   shot    into  the  cathode   surface  or   adsorbed  on  the  surface 
formed  by  the  sputtered  material.      Variation  of   sputtering  rate   is 
functionally  dependent  upon  a  negative  power  of   the  pressure,    as  both 
ion  velocity  and   sputtered  particle  free  paths   are   increased    in  this 
case.      3 

Acton   and  Swift   note  that   this   power   seems  to  be   approximately   -5, 
and  on  this  basis  they  derive      3 

dp     _     -d 
dt  p5 

P   =    (P06   -   6  ct)1/6 

A  tube  losing  5%  of  its  pressure  in  two  hours  will  lose  20%  of  its  pres- 
sure in  the  next  four  hours.   This  effect  is  known  as  cleanup  and  can  be 
employed  as  a  means  of  purifying  the  gas  in  a  discharge  tube.   All  ele- 
ments other  than  helium  have  ionization  energies  which  are  lower  than 
that  of  neon,  and  will  therefore  be  ionized  preferentially  in  any  region 
of  the  discharge.   The  ions  move  to  the  cathode  region  under  the  influence 
of  the  electric  field,  where  cleanup  takes  place,  leaving  the  body  gas 
free  of  impurities. 

Druyvesteyn  and  Penning  note  that  the  Faraday  dark  space  contains  a 
"red  aureole"  in  very  pure  neon.  1 16  J  This  is  believed  to  be  due  to  the 
excitation  of  metastable  atoms  which  have  diffused  into  this  region  from 
the  positive  column  or  the  negative  glow.   It  was  used  as  a  check  on 
purity  of  the  discharge  gas  in  all  experiments.   A  mass  spectrometer  gives 
the  only  certain  test  of  the  constituent  gases  of  the  discharge,  and  is 
impractical  for  ordinary  discharge  work. 
C.   Positive  Column 

Macroscopic  properties  of  the  plasma  of  the  positive  column  such  as 
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the  complex  conductivity  and  the  spectrum  of  the  emitted  light  are  deter- 
mined by  the  microscopic  processes  governing  thermal  equilibrium,  ioniza- 
tion, mobility  and  diffusion.   The  development  of  moving  striations  is 
dependent  upon  each  of  these  factors. 

In  this  section,  after  a  brief  description  of  plasma  and  the  Debye 
length,  thermal  equilibrium  is  discussed.   Elastic,  inelastic  and  electron- 
electron  interactions  are  included  since  each  of  these  is  important  in 
determining  the  electron  energy  distribution  function.   The  form  of  this 
function  determines  whether  ionization  will  be  predominantly  direct  or 
stepwise,  each  case  giving  rise  to  moving  striations  with  different 
characteristics. 

The  diffusion  of  charged  particles,  formation  of  space  charges  and 
the  effect  of  electric  field  changes  upon  the  temperature  of  electrons 
are  all  processes  which  have  been  studied  as  a  result  of  theoretical 
papers  on  moving  striations. 

1.   Plasma 

Before  proceeding  with  this   description  of   the  discharge,    it 
should   be  pointed  out  that   the  weakly   ionized   gas   of   the   positive  column 
is   a  plasma.      Irving  Langmuir   first   applied  this   term  to   the  positive 
columns   of  discharges   with  which  he  was   working   in  order  to   convey  the 
idea  that   positive   and  negative   charge  densities   were   approximately 
equal.      Poisson's    equation 

v   •    E   =   4TrqofN+   -   N_  ]  (esu)  (59) 

states  that  large  electric  fields  will  be  set  up  whenever  N+  /  N_,  since 

10   -  ? 
the  values  of  N+  are  of  the  order  of  10   cm  °.      The  mobility  of  the 

charged  particles  will  cause  this  imbalance  to  be  adjusted.   The  restor- 
ing motion  will  also  cause  oscillations,  but  characteristic  frequencies 
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where 


for  oscillations  of  both  ions  and  electrons  are  far  above  the  range  of 
this  experiment.   This  internal  struggle  for  neutrality  gives  rise  to  a 
characteristic  distance  for  the  decay  of  an  artificially  introduced 
potential  such  as  at  a  probe  or  the  tube  wall.   Known  as  the  Debye  length 
or  Debye  shielding  distance,  it  is  given  by 

\   =   ! (60) 

k  =  Boltzmann's  constant 
Te  =  electron  temperature,  degrees  Kelvin 
n   =  electron  density,  cm"J 
qQ  =  electronic  charge,  esu  . 

It  must  be  added  that  all  extraneous  fields  are  not  eliminated, 
and  the  plasma  does  not  succeed  in  eliminating  local  electric  fields  at 
all  times.   In  fact,  the  question  of  how  the  positive  column  supports 
a  macroscopic  polarization  such  as  that  which  accompanies  moving  stria- 
tions  is  the  very  question  which  has  perplexed  investigators  for  years. 
This  will  be  clarified  in  the  section  on  wave  of  stratification  theory. 

2.   Thermal  Equilibrium 

The  mechanisms  involved  in  attaining  thermal  equilibrium  are 
rather  complex,  since  the  electric  field  acts  as  an  energy  source  for 
the  system,  directly  affecting  both  ions  and  electrons,  while  these 
charged  particles  exchange  energy  among  themselves,  with  each  other  and 
with  the  neutral  gas.   Finally,  there  is  a  large  energy  sink  for  high- 
energy  electrons  in  the  form  of  inelastic  collisions  leading  to  excita- 
tion and  ionization.   Several  attempts  have  been  made  to  treat  this 
problem  theoretically,  and  as  each  work  adds  more  insight  into  these 
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fundamental  processes,  the  important  aspects  of  each  are  given. 

Since  the  maximum  fractional  energy  which  a  particle  of  mass  m~ 

can  lose  in  an  elastic  collision  with  a  particle  of  mass  m-,  is 

^m-ioin 
1  -  '  •  -  £   -  '.'  j  an  electron  will  lose  a  small  fraction  of  its  energy  in 
(m-,  +  riv?)^ 

colliding  elastically  with  a  neutral  neon  atom  -  approximately  1/37,000 
or  2.7  x  10"""  -  while  an  ion  will  lose  up  to  1/2  its  energy  in  colliding 
with  a  neutral.   This  allows  the  electrons  to  retain  more  energy  gained 
from  the  electric  field  than  the  ions.   Typical  average  electron  ener- 
gies are  of  the  order  of  two  electron  volts  (eV) ,  equivalent  to  a  kinetic 
temperature  of  about  23,000  degrees  Kelvin.   Ions  also  gain  energy  from 
the  electric  field,  but  ion  temperatures  differ  little  from  room  temper- 
ature due  to  the  efficient  transfer  of  energy  to  the  background  gas. 
Thermal  equilibrium  does  not  exist  between  these  three  species.   The 
interaction  between  them  and  among  particles  of  a  single  type  all  con- 
tribute to  establishing  a  steady  state  distribution  of  particles  with 
respect  to  energy.   The  average  energy  of  this  distribution  is  then, 
of  course,  the  energy  associated  with  the  temperature. 

The  distribution  of  electrons  with  respect  to  energy  changes  as 
electric  field,  elastic  collisions,  electron-electron  interactions, 
and  inelastic  collisions  with  neutrals  are  taken  into  consideration. 
In  the  field-free  case,  the  equilibrium  distribution  is  a  Maxwell- 
Boltzmann  function,  with  the  temperature  of  each  component  equal  to 
that  of  the  walls.   When  an  electric  field  is  added  with  only  elastic 
collisions  allowed,  and  the  Coulomb  interaction  between  charged  particles 
is  ignored,  the  distribution  function  is  known  as  a  Druyvestyn  distri- 
bution.  16   A  comparison  of  these  two  forms  of  the  distribution 
function  shows  that  the  effect  of  elastic  collisions  with  the  neutrals 
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f  (e)  ■  -Ar Ne  [*  e"x  dx  ]  (m-b)  «» 


*  <^    -  2  r    k     -0.55x^,   I 

£(€)    -   p(   3  )        Ne|^x2e  dxj  (62) 

(Druyvestyn) 
x   =     C/€q 
€     =   average  electron  energy 

is   to   reduce  the  number  of   high  energy  electrons.      These  high  energy 
electrons    are  essential  to   excitation  and   ionization  from  the  ground 
state.      Thus   the  Maxwell-Boltzmann  function  will   cause  more  excitations 
and   ionizations   for   a  given  temperature. 

Druyvestyn  also  derived  the  distribution  function  for  the   case   in 
which   inelastic  collisions   are   allowed,   and  Smit's   calculations   based 
on  this   work  showed  the  decrease   in  high  energy  electrons  which  can  be 
expected   in  helium.      73       This   loss   is   reduced  as   electric   field   is 
increased. 

Cahn  considered  elastic  collisions   and  charged   particle   interactions. 
8        His   results   showed   that   an  increase  of   electron  density  caused   a 
change  from  the  Druyvestyn  toward  the  Maxwellian  form  of  the  distribu- 
tion.     A  similar  result  was   obtained   by  Gurevich.   I  23  j 

Kagan  and   Lyagushchenko   appear  to  be  the  first    to   include  both 
elastic  and   inelastic   collisions   and  electron-electron   interactions   in 
a  calculation  of   the  distribution  function.      28        Their   result   for   a 
pressure  of  five   torr  and   200  mA  is   graphed  against   the  Maxwellian 
function  and  shows  that   the   number  of  particles    at   energies   above   9 
electron  volts    is   an  order  of  magnitude   smaller  than  the  number  for  a 
Maxwellian  distribution  corresponding  to   the  same   average  energy.    J  29  I 
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In  each  of  these  cases  it  is  evident  that  elastic  collisions  with 
neutral  particles  cause  the  deficiency  of  high  energy  particles.   Re- 
duction of  electric  field  and  electron  density  both  serve  to  make  the 
distribution  function  short  of  high  energy  electrons.   Thus  non-Maxwel- 
lian  electron  velocity  distributions  are  expected  for  pressures  and 
currents  in  the  experimental  range,  with  the  departure  from  the  Maxwal- 
lian  form  becoming  larger  as  pressure  is  increased. 

The  distribution  function  of  electrons  in  neon  has  been  measured 
by  Vorob'eva,  Kagan  and  Milenin  over  the  pressure  range  0.06  to  1.6 
torr.   82   At  0.06  torr  and  200  mA,  the  distribution  was  essentially 
Maxwellian.   For  a  pressure  of  1.0  torr  and  100  mA,  the  fraction  of 
electrons  with  energy  centered  at  1 7  eV  is  half  that  corresponding  to  a 
Maxwellian  function  .   This  shows  the  change  of  the  function  brought 
about  by  a  pressure  change.   Data  taken  at  the  same  pressure,  and  25  mA 
discharge  current,  show  that  the  distribution  decreases  more  rapidly 
with  energy  than  the  function  obtained  at  100  mA.   Electron-electron  i 
interactions  tend  to  make  the  distribution  function  Maxwellian.   These 
interactions  are  reduced  as  current  is  decreased,  and  this  explains  the 
further  departure  of  the  electron  distribution  function  from  Maxwallian. 

3.   Ionization  and  Recombination 

Ion  pairs  are  produced  in  the  positive  column  mainly  by  in- 
elastic collisions  of  electrons  with  neutral  atoms  which  are  either  in 

the  ground  state  or  in  a  meta-stable  state.   Photoionization  plays  a 

o 
minor  role,  as  ultraviolet  photons  with  wavelength  less  than  2400  A 

are  required  for  ionization  from  metastable  states,  and  vacuum  ultra- 
violet photons  are  required  for  ionization  from  the  ground  state. 

The  general  equation  for  the  rate  of  ionization  by  electrons  can 
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be  written 
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e  Ul 


oo 

(63) 


+   E  Ns    J        Q8iCu)-«f0(u)duj 

s  uD.  J 


The  first    integral   represents   direct   ionization   (from  the  ground   state) 
while  the  second   represents    ionization  from  excited   states.      Summation 
is   over  all  excited   states,    although  only  the  four  states  with  config- 
uration  2p^3s  have  populations   which  are   significant.      Two  of  these   are 
metastable.      Other  symbols   not   previously  used   are 

u  =     dimensionless   electron  energy      (eV/kT  ) 


ui 


=     eV./kT   :    V-    is  the   ionization  energy  of   neon 
Qj(u)      =      ionization  efficiency  for  electron  energy  u; 

approximated   by  N.Qj^   =   0.055  (V  -   20.6) 
Ng  =     density  of  excited   atoms   at  the   level  denoted 

by  s 

u  =     e(V.    -  Vs)/kT 

s-  i  &  e 

N«Qi   =   y.0.055(V  -  V-  -  Vs);  /is  determined  experi- 
s 

mentally,  corresponds  to  a  different  cross  section 
for  the  respective  states 
N      =  neutral  particle  density 
The  equation  for  Q.  expresses  the  fact  that  ionization  cross  section  of 
the  electrons  increases  with  increasing  electron  energy.   This  is  true 
for  the  energies  which  electrons  can  obtain  under  the  conditions  of 
these  discharges. 

In  view  of  the  deficiency  of  high  energy  electrons  noted  in  the 
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previous  section,  direct  ionization  would  be  expected  to  contribute 
fewer  ion  pairs  than  stepwise  ionization.   A  calculation  of  the  rela- 
tive importance  of  the  two  methods  of  ionization  shows  that  stepwise 
ionization  provides  more  ions  by  three  orders  of  magnitude  than  direct 
ionization.   29   Comparison  with  a  Maxwell -Bo It zmann  distribution  shows 
that  the  distribution  function  derived  with  electron-electron,  elastic, 
and  inelastic  collisions  taken  into  account  produces  fewer  ions  from 
direct  ionization  by  four  orders  of  magnitude.   These  results  lead  the 
authors  to  conclude  that  "the  role  of  direct  ionization  is  negligibly 
small  compared  with  that  of  stepwise  ionization".   30   This  is  a  con- 
clusion which  must  be  '^contrasted  with  experimental  findings  on  moving 
striation  behavior. 

The  change  of  light  emitted  by  the  discharge  can  be  considered  pro- 
portional to  the  perturbation  of  electron  temperature.   44   Equation 
(63)  shows  that  the  change  of  ionization  rate  can  also  be  considered 
directly  proportional  to  the  change  of  electron  temperature  for  small 
perturbations.   These  changes  of  light  emitted  and  ionization  rate  with 
electron  temperature  are  important  to  mathematical  formulation  of  the 
striation  processes,  and  to  interpretation  of  experimental  results. 

Recombination  of  electrons  and  ions  in  the  volume  of  the  positive 
column  is  negligible  compared  to  the  rate  at  the  tube  walls.   Consider- 
ing the  requirements  for  conservation  of  energy,  linear  momentum,  angular 
momentum  and  spin,  it  is  evident  that  the  surface  atoms  of  the  tube  wall 
provide  much  better  opportunities  for  these  conditions  to  be  satisfied. 
It  is  thought  that  ions  and  electrons  arrive  at  the  walls  separately  and 
the  recombination  takes  place  through  motion  along  the  walls.   37  I 
Since  this  seems  to  have  no  bearing  upon  the  processes  of  the  discharge 
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involved  in  these  experiments,  it  will  not  be  considered  further. 
4.   Mobility  and  Diffusion 

When  electrons  and  ions  move  through  a  neutral  gas  as  the  re- 
sult of  an  electric  field,  their  average  drift  velocities  are  given  by 

The  mobility  constant  for  ions  /  electrons,  IA  ,  is  a  function  of  the 
mean  free  path  and  average  thermal  velocity  40 

j j  e  A+  M+  +  M  t 

^  =  0-75l^-(-!r->2 

=  1.06  e  \   /   M+  c7  (64) 

^U     =  0.75  e  ).  /  mcT   .  (65) 

These  are  results  of  derivations  by  Langevin,  and  the  ionic  mobility 
is  quite  satisfactory.   Loeb  I  37   has  improved  the  electron  mobility 
expression,  but  the  formula  given  above  provides  values  which  are  satis- 
factory to  an  order  of  magnitude,  and  has  the  advantage  of  being  more 
easily  understood.   Using  values  of  mean  free  path  and  average  electron 
energy  given  in  chapter  IV,  the  following  values  of  mobility  can  be 
calculated: 

yO*+   =  1.24  x  104    (cm2 /volt  sec  ) 

4  9 

yU_   =   80  x  10       (cmVvolt  sec  ) 
Note  that  the  electron  mobility  is  much  higher  than  the  ion  mobility. 
This  is  due  primarily  to  the  lower  electronic  mass,  as  may  be  seen  from 
a  comparison  of  the  two  mobility  formulas. 

Mobility  is  also  important  in  the  diffusion  processes  of  the  posi- 
tive column.   Since  most  recombination  of  ions  and  electrons  produced 
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in  the   column  takes   place  at   the  walls   of  the  discharge  tube,    a  radial 
concentration  gradient  exists    in  the  tube.      Electrons   and   ions   tend   to 
diffuse   toward   the   walls    with  their  motion  governed   by   concentration 
gradients   and   electric   fields    according  to  the  equation 


|~^   =    -    V  (D+  N+)    +  N+ yu+  E    . 


(66) 


The  opposite  signs  arise  from  the  difference  in  signs  of  the  charges. 
In  the  steady  state,  particle  flux  of  each  type  must  be  equal,  and  the 
plasma  condition  forces  N  ps   N  .   If  the  electron  temperature  across 
the  tube  is  constant  as  is  essentially  the  case,  D+  is  no  longer  a 
function  of  position.   17   The  electric  field  can  then  be  eliminated 
from  the  equation  giving 

v+  =  £- i-ii.  •  ^- —  (67) 


/<*    */*. 


N 


=  "  Da  '  ^L  (68) 


Diffusion  of  this  type  is  called  ambipolar  diffusion,  and  if  T  >>  T- 
(as  ig.  the  case  in  this  work),  the  coefficient  of  ambipolar  diffusion 
can  be  simplified   17   to 

Da  -  -y-1  (69) 

Mobility  and  diffusion  play  an  important  role  in  forming  local  space 
charges  in  the  positive  column,  which  then  change  the  electric  field, 
electron  temperature  and  in  turn  the  ionization  rate.   Pekarek  has  in- 
vestigated this,  beginning  with  continuity  equations  for  ions  and 
electrons  and  the  Poisson  equation.   53   By  assuming  perturbations  of 
the  form 
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N±   =   No±      +      n±       ln'^   No  (70) 

E    =   E0    +  e  |e|^   EQ  (71) 

these  equations  can  then  be  written 

(}n+  A^n+  An.  ^p 

ST  '^-fe-Wiir  -**>■&-*'  (72) 

dn-                   A2n                           <^n-  Ae 

— v =  D     — — h-  +  Z4     Eo  —5 +    /<    N~  t-  +  Q  (73) 

£t  ~     ^x2      /-    °    <3  x        A*-    0-57-      v        UJJ 

^e 

0   =  ^— -  -  4irqo-(n+   -  n_).  (74) 

Since  ion  concentration  changes  much  more  slowly  then  electron  concentra- 
tion, the  ion  concentration  is  assumed  given.  Substituting  equation  (74) 
into    (73)    results    in  the  differential   equation 

-*T  ~  D-  ~T^- +  ^"E°  "ST"  +^-V^on- 

=  Q   +  4TTQ0  NQ   n+  <75> 

Under  the  assumption  that  Q  =  0  and  a  Dirac  £ -function  initial  ion  per- 
turbation, it  is  shown  that  the  relaxation  time  for  the  electron  distri- 
bution to  be  established    is   given  by 

D_  (76) 

For  the  conditions  of  this   experiment 

X^       =        5.0  x   10~5   cm 
D_      =        1.8   x   10 6   cm2/sec 

2T_     =        2.8  x   10"11   sec 
Characteristic   frequencies    for  moving   striations   are  of  the  order  of 
0.5   -   50   kHz.      Electron   relaxation  time   is   several   orders  of  magnitude 


73 


shorter  than  times   associated  with  changes   of   striation  parameters.      The 
electron  density  thus    follows   that   of   the    ions    rapidly  enough  to   consider 
it   quasi-static.      With  this    result,    the  electron  density  can  be   written 
in  terms  of  the   ion  density   as 


where 


,00 

n_(x,    t)    =    /      G_(x  -|)    •    n+    (J,    t) 

-L  co 


(77) 


and 


G    (x)    =-4-«    e"X/£l    .    e-lX'^D  (78) 

2*g 


Ap        =       Debye   length      (eq.60) 
\        =        2D-^-'Eo 


(79) 


'  l  1     h 

*V        =         (-TT-     +     -JT>  (80) 

*D  1 

and   the   space  charge   is  then  given  by  the  equation 

Ax,    t)    =  qQ    J[S(  x    -  J   )    -   G_(x  -J)]    -    n+(|,    t)dj    .    (81) 
Equation    (81)    shows  that   the   space   charge   is  different   from  zero  when- 
ever  n+  ^  0.      Its  magnitude  depends   upon  n+   and  the  plasma  parameters 
through  equation    (78)  • 

There   is    another  consequence  of   the   ion  perturbation.      The  space 
charges   which   are  formed   create  electric   fields   which  perturb   the  main 
dc   field.      Since   the  electric   field  determines   the  electron  temperature 
through  the  equation 

~6\        3-       lfi     u2 


e      _      3  ^  16       Ue        R* 


2        •    —     E.>2  (82) 


where 


e 


U  =  electron  energy   in  eV 

^  =  electronic  mean  free  path 
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•k 

K    =   mean  energy  loss  of  the  electron  in  an 
elastic  collision  , 
the  change  of  electric  field  also  changes  the  electron  temperature. 
Pekarek  analyzed  this  change  as  a  function  of  distance  by  introducing 
perturbations  of  electric  field  and  electron  temperature  through  the 
relations   51 


U   =  Ue   +  6     9  «Ue  (83) 

e      eQ        i  i    eQ 

E   =  EQ  +  e     |  e| «  EQ  (71) 


o 
The  spatial  change  of  electron  temperature  becomes 


de    _ 


"3 


X 


=  axe  -  bxe  (84) 


al 

3TT 

EQ^ 

to  K 

bl  = 

'  3 
2 

16 

+  . 

37T 

"e"     1  a 

*o\  J    ^ 

if  second  order  terms  are  neglected.   In  arriving  at  this  relation,  the 

quantity  Ue   is  equated  to  the  energy  gained  from  the  electric  field 
o 

in  a  mean  free  path.   The  constants  a-,  and  b..  are 

-!£-  Ue°„   (2K*  +  UP.   <^K*  )  (85) 

(86) 

The  importance  of  equation  (8^)  is  the  shift  of  electron  temperature 
maximum  from  electric  field  maximum.   Pekarek  has  shown  that  this  spa- 
tial phase  difference  is  essential  to  amplification  of  the  wave  of 
stratification.   51 

Since  a  change  of  electron  temperature  also  changes  the  light 
emitted  by  the  discharge,  variations  in  light  emitted  along  the  dis- 
charge can  be  attributed  to  temperature  gradients,  local  changes  in 
electric  field,  and  ultimately  to  regions  of  local  space  charge. 
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The  formation  of  space  charges  and  electron  temperature  gradients 
also  changes  the  diffusion  processes.   During  the  derivation  of  the 
ambipolar  diffusion  coefficient,  it  was  assumed  that  the  electron  dif- 
fusion coefficient  was  constant  as  electron  temperature  is  nearly  con- 
stant across  the  tube.   A  temperature  gradient  along  the  column  means 
that  a  perturbation  of  ion  density  will  decay  through  electron  tempera- 
ture diffusion  as  well  as  through  ambipolar  diffusion.   In  order  to 
compare  the  magnitudes  of  these  terms,  Pekarek  and  Krejci   55   write 
equation  (66)  in  its  full  form 

I-     =  (D_)  VN_  +  (N_)  VD_  +  E^w_N_ 
and  denote  the  ratio  of  the  terms  of  interest  by 

R.  N-(VD-}  (87) 

D_(  VN_) 

noting  that   the  diffusion  previously  discussed  occurs   when  R-<^  1.      Ex- 
panding the  diffusion  coefficient    about  the  equilibrium  temperature 
and  making  use  of  the  relations 

D_   =   kT.yW.   /qQ  (88) 

6      =  k(T_   -  T_Q)  (89 ) 

the  expression   for  R  becomes 

r.  N-°ye    y  (90) 

kT_0  7n. 


where 

/4_0        dkT_        T_c 


(91) 


Except   for  the   case  of  elastic  collision  cross-section   increasing 
linearly  with   electron  energy,    (making  7-  0),    -/is   of  the  order  of  one, 
and  may  be   positive  or   negative.      Thus   determination  of  R  depends    upon 
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the  quotient  VQ/  7n.»  Equation  (77)  gives  the  form  of  the  electron  den- 
sity curve,  and  a  similar  equation  describes  the  electron  temperature  as 
a   function  of  position 


Gfl    (x  -  § )    n+    ( f )   d  § 


where 


/oo 
Ge 
-  oo 

G0(x)P=  G°   =   2  7rq0-b1    .   |__£_     e"X'  VU(x) 


(92) 


/7t>    rrV        N  2al     ^D  anX         ,  ' 

D.U(-x)    -  ; =-  e    l    -U(-x) 

IT)  l-(a^D)2 


*  D 

+    — T 5 e 

1-a 


(93) 


Calculating  the  gradients   of  6,   n.for  a  sinusoidal    ion  density  perturba- 


tion with  wavelength  much   larger  than  X.  gives   the  result 


R  = 


!<•■/, 


(94) 


An  earlier  paper  showed   51 

b4  =  3qQ/2  . 
The  value  of  R  is  thus  of  order  one.   This  means  that  diffusion  of 
electrons  by  temperature  gradient  is  of  the  same  order  of  magnitude  as 
diffusion  due  to  a  concentration  gradient. 

Investigation  of  these  discharge  processes  has  come  about  as  a  re- 
sult of  research  into  moving  striations.   Some  are  applied  to  this  work 
to  explain  experimental  results.   Others,  such  as  the  diffusion  of 
electrons  due  to  a  temperature  gradient,  have  not  yet  been  applied  to 
striation  work. 
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APPENDIX  II 
WAVES  OF  STRATIFICATION 

This  section  gives  a  brief  description  of  waves  of  stratification 
and  the  various  theories  which  have  been  set  forth  before  reporting  work 
on  irradiation  of  the  discharge  by  an  external  source.   Parameters  of 
acoustically  excited  waves  of  stratification  are  also  presented. 
A.   Theory  of  the  Wave  of  Stratification 

Several  attempts  to  describe  mathematically  the  development  of  the 
wave  of  stratification  on  the  basis  of  the  microphysical  processes  of 
the  positive  column  have  been  set  forth.   The  earliest  of  these  works 
was  published  by  Pekarek  in  1957.   *+3   In  this  approach,  the  observed 
structure  of  the  wave  is  postulated,  the  changes  of  electric  field  and 
space  charge  are  related  to  this  structure,  and  the  production  and  decay 
of  the  space  charges  are  included  by  relaxation  times  characteristic  of 
the  processes  involved.   Each  maximum  and  minimum  of  light  intensity  is 
indexed  so  that  the  light  area  numbered  2k  is  actually  the  k   light 
area  (figure  1-1).   Perturbations  of  steady  state  quantities  are  denoted 

by 

9   electron  temperature 

n   charged  particle  density 

m   met as table  atom  density 

q   space  charge  density, 
and  the  constants  a<  ,  /&   ,  are  introduced  to  describe  the  decay  of  space 
charge  to  zero  and  the  removal  of  the  metastable  perturbation  by  ioniza- 
tion. 

With  these  quantities  describing  the  column,  and  assuming  the  light 
intensity  changes  to  follow  those  of  the  electron  temperature,  the 
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periodic  changes  in  the  plasma  are  described  by  two  chains,  depending  on 
whether  the  ionization  is  predominantly  stepwise  or  direct.   Both  chains 
begin  with  a  change  of  electron  temperature  which  is  assumed  to  be  ac- 
companied by  a  change  in  ionization  rate: 

I.  +n m  +q  *>   0  -0. 

S       ^ ' 

II.  +m ^  +n      -  +q ^   0        -0 

In  chain  I,  direct  ionization  is  assumed  to  produce  an  increase  in  the 
number  of  charged  particles  of  both  signs,  and  space  charges  form  due  to 
the  higher  mobility  of  the  electrons.   These  space  charges  decay  to  zero 
with  relaxation  time  tl    =  \/o<,   but  cause  a  change  in  electron  tempera- 
ture of  opposite  sign  in  the  adjacent  region  on  the  anode  side.   In 
chain  II,  metastable  atoms  produced  by  the  increase  in  electron  tempera- 
ture are  ionized  in  a  characteristic  time  7^,  =  l/£  •   From  this  stage, 
the  processes  proceed  as  in  chain  I. 

It  is  convenient  to  use  a  coordinate  system  in  which  the  striations 
are  stationary  for  writing  the  equations  which  describe  these  processes, 
and  transform  the  solution  back  to  the  laboratory  frame  (in  which  the 
striations  are  moving)  after  the  solution  is  obtained.   Using  this  method, 
the  following  equations  can  be  written: 

Chain  I 


dn 


8-   =  A'e,   -  <*ny  (95) 


Chain  II 


dt 


i*f_=  B*9y   -/9  m,  (96) 


dn 


id_  =  i  my   -  ^n^  (97) 
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For  both  chains,    the  electron  temperature   is   assumed   to   follow  any  change 
of   electric   field    immediately.      The   equation  expressing  this    is 


8y(t)    =   "  |V/    (t)   =   -k''    Vl    (t) 


(98) 


The  second  equality  states  that  the  space  charges  are  formed  whenever 
there  is  a  change  in  charged  particle  density.   Substituting  this  equa- 
tion into  the  equations  describing  chains  I  and  II  gives  the  following 
differential  equations: 

d6. 


V 


+  o<e^   =  -ckA0  (99) 


dt     '      y-i 


d2ey 


+  (  o<  +  Q   )  JSlL.  +  *(&,      =  -Bo^_!  '     (100) 
dt 


These  equations  are  solved  under  the  initial  condition  of  a  Dirac  delta 
function  increase  of  the  electron  temperature,  and  the  solution  gives  the 
time  development  of  electron  temperature  in  each  of  the  numbered  regions. 
The  time  of  occurrence  of  the  maximum  amplitude  of  each  region  is  deter- 
mined, as  are  the  times  of  occurrence  of  half  maximum  intensity.   Since 
the  time  of  occurrence  of  the  maximum  of  light  intensity  is  related  to 
the  propagation  velocity  of  the  wave  ( u)  and  the  striation  velocity  (v)  , 
both  of  which  are  measurable  quantities,  and  are  related  in  turn  to  the 
relaxation  times  through  the  striation  spacing  or  wavelength  (  )0  »  the 
relaxation  times  can  be  checked  from  the  experimentally  determined  values 
for  agreement  with  the  hypotheses  as  to  the  governing  processes  in  each 
chain.   This  determination  of  relaxation  time  has  been  a  most  useful 
result.   It  can  be  shown  that 


T°  "  ~ uTT"  (101) 


80 


where  the  relaxation  time  is  related  to  the  plasma  processes  through 

C°        "ST  +  -J-  1     2  (102) 

for  cases    in  which  there  are   two   processes.      Also  predicted  by  this 
theory  are 

(a)  oddness   of  the  wave  of   stratification   (initial   conditions  of 
opposite   sign  produce  waves   of   stratification  with  light   and  dark  areas 
reversed) 

(b)  a  spatial   decrease   in   amplitude  of   the  wave  group  followed   by 
an   increase  when  the  amplification   is    larger  than   1.0 

(c)  direct   dependence  of  wave  of  stratification  amplitude   upon  the 
amplitude  of  the   initial  disturbance 

(d)  growth  of  width  of   the  packet    in  time  according  to  t 2. 
These  predictions  were  verified   in  a   later  experimental  work,   with  the 
exception  of  the  width  of  the  packet   of  the  fast  wave.      48        The  experi- 
mental results    in  this   case   agreed   with  theoretical   predictions    in  form 
although  not    in  magnitude. 

It   should  be  emphasized   that   this   theory  takes   the   layered  structure 
as    a  basis,    attributes  the   light    intensity  variations   to  changes    in 
electron  temperature,   and   attributes   ensuing  changes    in  electron  tempera- 
ture to  the   space  charges   produced  by  the  higher  mobility  of  the  electrons. 
Basic  microphysical  processes    in  the  positive   column  are  not   used   to 
predict   the   response  of  the  plasma  to  the  perturbation  of    ionization 
rate. 

A  later   theoretical  work  uses    a  simplified   form  of  Poisson's   equa- 
tion and  the  equation  of  continuity   for   ions   with  only  production  of   ions 
due  to   changes    in  the   ionization   rate  being  considered.    1 51  |     The   Debye 


81 


length  is  essentially  assumed  infinite,  as  the  electrons  arising  from  the 
perturbation  of  the  ionization  rate  are  assumed  to  diffuse  immediately 
upon  being  created.   Perturbation  relations  of  the  form 

E  =  EQ  +  e        |e|  «E  (71) 

N  =  NQ  +  n        |n|  «N  (70) 

substituted    into   these  equations  give  the  equations 

■g^-  =   4TrVV  (esu)  (103) 

A^-  =   2'Nne  (104) 

6t 

where  g      is   the  rate  of   change  of   ionization  rate  with  electric   field, 
and  qQ    is   the  electronic  charge.      The  solution  to   this   system  of   equa- 
tions under  the   initial   condition  of   a  unit   step   in   ion  density   in  both 
time  and  space   is   a  Bessel   function  with  the  product    (xt) 2  as   the  argu- 
ment.     With  coordinate  axes   chosen   such   that   positive  X  is    in  the  direc- 
tion of  the  electric   field,    and  axes   designated   as    in  the  time-space 
diagram  of  figure    1-1    (negative  X  in  the  vertical   direction),    periodic 
solutions   arise   in  the  form  of   a  family  of  hyperbolae    (figure  ik-2)  .      If 
very   short   and  very   long  wavelength  striations   are  excluded   from  con- 
sideration by  drawing  the  dotted   lines    shown  approximately   15  degrees 
either   side   of  the  quadrant   bisector,   the  resulting    lines   have  the  gen- 
eral  appearance  of   the  wave   of   stratification.      This  diagram  then  shows 
the   propagation   in  the  cathode  to   anode  direction,    the  motion  of   stria- 
tions   in  the  cathode  direction,    and  convex  dispersion  of  the  wave    (dis- 
cussed   later) . 

This   simple   set   of   equations  thus   contains   the  necessary  processes 
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for  the  propagation  of  a  periodic  structure  with  the  properties  of  the 
wave  of  stratification,  provided  short  and  long  wavelengths  are  eliminated, 
An  improved  treatment  must  include  these  processes.   Pekarek  lists  the 
three  important  physical  processes  which  cause  the  wave  of  stratification 
as 

(a)  dependence  of  the  idnization  coefficient  on  the  electric  field 
in  the  plasma 

(b)  the  production  of  space  charges  due  to  different  diffusion 
rates  of  the  charged  particles 

(c)  creation  of  additional  electric  fields  due  to  the  creation  of 
space  charges. 

Since  perturbation  of  the  electric  field  causes  a  perturbation  of  the 
electron  temperature,  the  ionization  rate  is  dependent  upon  electron 
energy.   Later  formulations  have  used  this  dependence  instead  of  the  de- 
pendence upon  electric  field.   Moreover,  the  shift  of  electron  tempera- 
ture maximum  from  the  electric  field  maximum  has  been  shown  to  be  essen- 
tial to  the  amplification  of  the  wave  of  stratification.  52  I  If  the  wave 
is  rapidly  attenuated,  its  characteristics  cannot  be  measured,  while 
large  amplification  causes  nonlinear  behavior  of  the  wave,  invalidating 
the  assumptions  of  the  linearized  theory.   Thus  the  behavior  of  the 
electron  temperature  is  important. 

Pekarek  and  Krejci  base  an  extensive  theoretical  work  on  the  equa- 
tions of  continuity  of  ions  and  electrons,  the  equation  for  electron 
temperature  variation  in  a  spatially  varying  electric  field,  and  Poisson's 
equation.   54,  56   These  are  written  as 
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Figure  4-2.   Solution  to  simplified  system  of  equations 

describing  electric  field  variations  in  the  positive 
column  following  unit  step  initial  conditions. 
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^r-   (n+   -   n_)    +   2q^No0  (105) 


— =  D       x    r>      +   A*    ^cPr—  +  /*  Eo     a? 

^t  3  z2  ~~    °<3z       ^"  °    <3Z 

+   3^   No0  (106) 

0   =  4^-  -  ^TTq    (n+   -  n  )  (74) 

di  ° 

0  =  -jg-  -   ai9  +  bxe  (84) 

Perturbations  of  the  form  given   in  equations    (70),    (71),    and    (83)    have 
been  applied.      The  terms   z  NQ9   and  -^ — (n+  -  n_)    represent   production 
of  particles  due  to   a  change    in   ionization  rate   and   loss  due  to   ambi- 
polar  diffusion.      The   term  z^No0  was  present    in  the   simple   set  of 
equation    (eqs.    0-02),    (103))    which  produced  the  time-space  periodic   struc- 
ture of  figure   4-2.      This    set    included   no  damping  mechanism.      Ruzicka 
and   Pekarek  carried  out   a   solution  to  the  set  of  equations    given  above 
with  the  term  z  N  9  deleted.      70        The   system   in  this   case  gave  only  a 


monotonic  decay  of  an  aperiodic  disturbance  <in   ion  density,    indicating 
that   the  change  of   ionization  rate  with  electron  temperature    is    a  neces- 
sary condition   for  the  excitation  of   a  wave  of   stratification. 

i  m     '    The  equations  are  simplified   for   solution   in  the  following 

manner: 

62n+ 
(1)  D+  —  is  neglected  as  the  changes  of  ion  density  due  to 

<3  z2 

additional    electric  fields  will   be   larger  than  the  diffusion  due  to   con- 
centration gradient 
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<3n       » 

(2)  ~   and  ZqNo0  are  neglected  in  the  electron  balance  equation 
as  the  change  of  electron  density  has  been  shown  to  adjust  to  the  change 
of  ion  density  several  orders  of  magnitude  faster  than  the  average  life- 
time of  the  charged  particles.  ["5  3  j  Electron  density  is  thus  quasi- 
stationary. 

(3)  the  system  is  transformed  to  a  coordinate  system  at  rest  with 
respect  to  the  ions  by  Z  =  x  +  ^  E  t 

(4)  the  term  /t/+  is  neglected  in  the  sum  /*  +  u_   . 
The  equations  which  result  are: 


-j£-  =  V+No  "jf-  " 


±r-   Cn+  -  n_)  +  z'nq9 


n  -  ™   c3  n-        „   O n-        „  (3  e 

0  -  D.  -j-5-  ♦  ^.E0  —  ♦  /-.><„—  a08, 

c*e 
0  =  -g^-  -  ^TTq0(n+  -  n_)  (74) 

0  =  -gJL  -  aie  +  bxe  .  (84) 

These  equations  receive  no  further  simplification  before  a  general 
solution  is  obtained.   The  processes  included  in  these  equations  are 

(a)  change  of  ionization  rate  with  electron  temperature 

(b)  creation  of  space  charges  due  to  different  free  diffusion 
rates  of  ions  and  electrons 

(c)  creation  of  additional  electric  fields  due  to  creation  of 
space  charges. 

These  were  shown  to  be  sufficient  for  the  formation  of  a  periodic  time- 
space  structure.   Also  included  are  terms  for 

(d)  the  displacement  of  electron  temperature  with  respect  to 
electric  field  (necessary  for  amplification  of  the  wave  of  stratification) 
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(e)  damping  of  the  striations   by  ambipolar  diffusion 

(f)  finite  diffusion  coefficient   of  electrons 

(g)  finite  mobility  of  electrons 
(h)    ionic  mobility 

(i)    effect   of  electron  density  on   ionization  rate. 
The  time  and   space  development   of  the   ion  density   in  terms  of  these 
basic   processes    is   obtained  by  using  the  2-sided  Laplace  transform 
described    in      61 1.      The  equation  which  results    is: 

■oo 

K( 
'-oo  ^ 


2l!±-(x,   t)   =/    K(x  -j).'n+(f,    t)-df  -  — L_  n+   .  (109) 


In  this  equation,  the  term  K(x  -J  )  is  a  sum  of  three  exponential 
functions.   The  equation  in  this  form  is  very  difficult  to  solve.   By 
assuming  that  the  Debye  length  is  nearly  zero,  the  authors  are  able  to 
make  simplifications  and  finally  arrive  at  the  following  equation: 

An,  kT_       32n+       ,      n 


/OO 
eai(z  -S>  u(_z  +|}  n+  (^  t)  d| 


oo 

e3i 
-  oo 


Their  analysis  of  each  term  on  the  right  hand  side  of  the  equation  is 
important  and  can  be  summarized  as  follows: 

(a)  the  first  term  is  a  diffusion  term  which  damps  the  very 
short  wavelength  striations 

(b)  the  second  term  makes  amplification  possible  (note  that 
b,  =  0  means  full  damping) 

(c)  the  integral  gives  an  oscillatory  character  to  the  solu- 
tion and  also  damps  the  long  wavelength  striations  (through  the  term 
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a^) ,    therefore 

(d)   there   are  no   redundant   terms   in  this   equation.      All   the 
major  physical   characteristics    are  present,    and   any   additional   terms 
can  only  modify  the   solution,   but   not   alter   its   basic   form. 

Solution  of   this  equation  under   a  Dirac  &  -function  disturbance  of 
ion  density   at   time  t   =  0   has   the  form 

n+(Z,   t)   =  -^-    f  exp  fi(kz  -a»(k)t)    +  0(k)tl   dk  (111) 


where 


-  oo 


k    =    2  TT/^ 

}       is   striation  wavelength 

A-iQk 
"><*>    =        2  2  (112) 

kl    +   al 

D_ 
A10   =   zebl(E°   +   al  -jj—  >  (113) 


0(k)    =   -k2D„   +   z'b,    -S= 

a  8     1      JA 


A10al 


6    l    /"-  k2    +   a2 


(114) 


Equation    (111)    is    a  Fourier   integral   over   all  wave   numbers.      The  term 
#(k)    is   present   as   an  amplification   (or  damping)   term,    and  thus   shows 
how  the  microphysical   processes    affect   striation  wavelength.      Amplifica- 
tion will  be  maximum  for 


dg(k) 

6* 


k=k1 


=  01   =  0    , 


and  this  gives 


where 


kx  =  +  ax(  *  -  1)^  (115) 


1-TT 


al  Da  M 


n  =  +  ^^  <«« 
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Real  values  of  k^  exist  for  0<)\<1,  and  when  used  in  equation  (1]2)  give 

A10ki  Da   a?  -i  i 

uj.  =  — 1U  l       =       at   T     (-L.  -  iys  (li?) 

1     kj  +  a2         n      i 

Wojaczek  has  discussed  a  method  of  solving  an  integral  of  the  type 
shown  in  equation  (111).  I  88  I  Functions  W(k)  and  0(k)  are  expanded  in 
a  Maclaurin  series  about  the  resonant  wave  number  (k,)  retaining  only 
first  order  terms.   This  is  valid  only  for  large  distances  from  the  point 
of  origin  of  the  wave  (large  Z)  ,  and  assumes  that  there  is  large  damping 
for  all  other  wavelengths.   Using  this  technique,  the  integral  in  equa- 
tion (111)  can  be  evaluated,  giving 

n_  f    (z-"it)2       1 

n+(Z,  t)  =: 2 exp — -i +  ^t 

(41T|B|  t)%      L      ^l*         x  J 

r                         (z  -  wxt)2  tit  n 

*  cos   lklZ  "  wlfc  +  — T^T "  ~f  J  (118) 


where 


„2  ,,2  „2  m2 

^1      +  0X  ^1      +  01 


b.    =  1 i —  q      = 

1 


1  ,,2  »2  %  _1      -6*f 

B      =  _L_   (cu'      +   j?      )2  =   tan      (— -^-  ) 

2  *  1  (^ 

The  authors  plot  the  solution  as  a  function  of  time  for  a  given 
coordinate,  noting  that  the  plot  has  the  typical  form  of  the  wave  of 
stratification  (figure  2-3(a)) .  Group  and  phase  velocities  obtained 
from  the  customary  relations 

v   =  ^L-     =   D     Jl-  <119> 

fcj        a   n 

ak         a  n        l  (120) 
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are  of  opposite  sign.      The  wave   has   anode-directed  group  velocity,    (due 
to  the  negative   sign)    and   is   thus   of   the  general   nature  of   the  wave   of 
stratification.      By   using  striation  wavelength  and   frequency  obtained 
experimentally,   together  with   an  estimated   relaxation   length    (1/a,)    for 
electron  temperature,    the  quantities   calculated   from  equations    (118,    119) 
and    (    120  )    are  found  to  agree  with  the  observed  values   of  these  to  with- 
in an  order  of  magnitude.      This   theory  therefore  provides    a  satisfactory 
description  of  the  wave  of  stratification  for   pressures   of   a  few  torr 
and  currents  of  order  of   an  ampere. 

Since  the  phase   and  group  velocities   are  of  the   same  order  of  magni- 
tude,   and    |u|^|v|,   this   theory   is    inapplicable  to  the  slow  and  fast 
waves   observed  at    lower  currents.      It   also   predicts   a  monotonic  decay  of 
light    intensity  on  the   cathode   side  of  the   initial   disturbance. 

Lee,   Bletzinger   and  Garscadden  have  solved  equation    (HO)   by  using 
a  Fourier  series   technique.      36  1    Their   results    show  one  significant 
difference,    in  that   a  forward  wave,    propagating   from  anode  to   cathode, 
is   predicted  whenever  ^^27T/a^.      This   has   recently  been  observed    in  ex- 
periments  on  a   35   torr   argon  discharge,    j  20 

Wojaczek  has  made  theoretical   and  experimental   investigations  of 
small   amplitude  moving   striations   and  waves   of   stratification,   j  84-86, 
88-90  I.      Using  conservation  equations   for  electrons   andiions,   Poisson's 
equation  and  the   equation  for   energy  balance  of  the   electrons    (including 
thermal  conductivity),   he  assumes    small   perturbations   of  the  variables 
with  plane  wave  time   and  space   dependence.      A  dispersion  relation   is  ob- 
tained,  and  with  the  assumption  of  a  complex  wave  number,    amplification 
as    a   function  of   frequency   is   derived.  I  88        Comparison  with  experiment 
shows   good   agreement   in  form  for  both  the   amplification  and  dispersion. 
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A  plot   of   R   •  U>Q   against   R   .    pQ,    where 
R        =     tube  radius 
cOQ     =      resonant   frequency 
pQ      =     pressure 
gives   a  curve  which   is    in  good   agreement   with  experiment   for  small 
values  of  R   •    pQ.      Pupp's   earlier  experiments   also   gave  data  which  fit 
this   curve.      66 

Dispersion   is  the  variation  of   phase  velocity  with  wavelength    (or 
frequency).      A  dispersion  curve   is   usually   constructed  with  the  fre- 
quency   (**>)    as   a   function  of  the  wave  number    (k   =   2  Tf/X) .      Figure   4-3 
shows   such   a  curve.      Since  the   phase  velocity   and  group  velocity  can  be 
calculated   from  the  relations 

v  =  cu/k  (119) 

u  =    dcu/dk,  (120) 

the  dispersion  curve   is   useful    in  determining  the  propagation  character- 
istics  of  the  wave  of   stratification.      Krejci  and  Pekarek  have  analyzed 
the   combinations  of  backward   and   forward  waves   and  their  associated  dis- 
persion curves.      34        In  their  work   it    is   shown  that   the  curvature  of  the 
striation  paths    also   determines   the  curvature  of  the  dispersion  curve  at 
the  point   of  the  resonant  wavelength.      Figure   44  shows    the    idealized   forms 
of  the  waves  of   stratification   in    (a)    neon,    and   (b)    a  neon-hydrogen 
("***  77l)   mixture.      The  ratio  of  group  and   phase  velocities    in  each  wave 
are   approximately   the  same;    magnitudes    are  different'. ~    .      Dispersion  is 
different   in  that   the   lines    are  curved   away  from  the  origin    (convex)    in 
(a),    and  toward  the  origin    (concave)    in    (b) .      The  portion  of  the  disper- 
sion curve  corresponding  to   each  of   these  cases   is    shown   in  figure   4-4. 
This  dispersion  curve    is   predicted  by  the   solution  to   equation    (111).      In 
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k  [cm"'] 


Figure  4-3.   Dispersion  curves  given  by 

(a)  simplified  model  with  infinite  Debye  length 

(b )  finite  Debye  length  theories. 
Diagram  after  Pekarek.  [6®} 
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Figure   4-4.      Time-space  development    of    light    intensity    in 
discharges    in    (a)    neon,    (b)    neon-hydrogen,    following 
an    impulsive    perturbation  of    ion  density   at    the    point 
X  =   0,    t    =   0. 
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most  cases,  however,  the  resonant  wavelength  is  such  that  the  convex 
portion  of  the  dispersion  curve  is  observed.   Curve  (a)  of  figure  4-3 
is  predicted  by  the  solution  to  equations  (103)  and  (104).   The  assumption 
in  the  latter  case  was  an  infinite  Debye  length,  while  in  the  curve 
designated  (b) ,  a  zero  Debye  length  was  assumed.   Thus  the  plasma  pro- 
cesses have  not  been  formulated  well  enough  to  allow  information  to  be 
derived  from  the  dispersion  curve.   The  authors  outline  a  method  for 
obtaining  data  from  the  photographs  and  constructing  dispersion  curves, 
and  this  both  sets  a  standard  and  gives  a  simple  method  of  accomplishing 
this  work. 

Pekarek  offers  an  explanation  of  the  concave  dispersion  in  the  neon- 
hydrogen  mixture,  and  of  the  anode-directed  striation  motion  in  pure 
hydrogen.   59,  60   The  basic  differential  equation  (eq.  110)  has  a  quad- 
ratic ionization  term,  and  in  the  case  of  a  finite  Debye  length,  this 
term  may  shift  the  ionization  maximum  toward  the  anode  from  the  point 
of  maximum  electron  temperature.   Since  the  ionization  maximum  is  shifted 
toward  the  anode,  electric  field  will  also  be  increased  in  this  direction, 
causing  the  temperature  to  increase,  and  the  chain  of  events  to  recur. 
This  is  the  tentative  explanation  which  he  offers  for  anode-directed 
striation  motion  in  hydrogen.   For  the  case  of  the  mixture,  contribution 
from  the  quadratic  ionization  term  increases  as  striation  amplitude  in- 
creases, giving  a  larger  anode-directed  component  to  the  striation 
velocity  at  the  larger  amplitude  phase  of  the  striation.   Thus,  the  case 
of  concave  dispersion  is  thought  to  be  associated  with  non-zero  Debye 
length. 
B.   Striation  Boundaries 

An  empirical  relation  for  the  upper  and  lower  striation  boundaries 
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is  given  by  Rutscher  and  Wojaczek.   69   Using  the  formulas  given  for  the 


upper 


and  lower 


Al 

p 


As 


Ls  =   p  -  Bs 

boundary  currents,  it  is  possible  to  make  calculations  for  the  tubes 
used  in  this  experiment.  (Subscripts  1  and  s  come  from  the  German  words 
"langsam"  and  "schnelle",  meaning  slow  and  fast,  respectively.  Slow  and 
fast  striations  usually  arise  at  the  respective  boundaries.)  Quantities 
are  defined  as  follows  for  neon,  with  p  given  in  torr  and  tube  radius  in 
centimeters: 

Radius 


as 
A„  = 


S       e< 

R 


Al 

Bl 

X 

8.0 

.1 

2.43 

6.6 

.07 

1.83 

6.4 

.06 

1.59 

as   = 

>    .24 

<*   =  .63 


s_ 

f3  Ds 


Bs  =  — IT  bc  =  .93 


(3   =  1.1 
For  tube  I,  R  =  1.0,  and  for  tube  II,  R  =  2.5.   Tube  I  was  used  in 
the  irradiation  experiments,  tube  II  in  the  acoustic  experiments.   If  we 
assume  the  following  values  for  tube  II,  which  can  be  obtained  from  this 
table  by  interpolation,  then  the  striation  boundaries  can  be  calculated 

A.  =6.5        ~X    =    1.71 
Bx  =  0.65 
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Discharge 

Current 

(mA) 


Neon   pressure    -   torr 

Figure   4-5.      Boundary   currents    for   striat ion-free 
operation   as    a   function  of   pressure 
in   tube   I.      Diameter:    1   cm,    length    30   cm. 
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for  each  of  the  tubes.   Results  are  plotted  in  figures  2-3  and  4-5  as 
dotted  lines,  and  the  experimentally  determined  values  are  shown  as  points 
with  a  curve  drawn  through  to  allow  comparison  with  the  results  of  this 
calculation.   It  is  evident  that  the  striation  boundaries  found  in  this 
experiment  are  much  different  from  those  predicted  by  the  empirical  rela- 
tionships.  In  all  cases  the  experimental  boundaries  occur  at  larger  cur- 
rents, and  are  as  much  as  an  order  of  magnitude  larger  in  the  case  of  the 
lower  boundary.   With  the  exception  of  the  one  cm  tube  at  the  upper 
boundary,  all  experimental  points  fall  in  such  a  manner  that  a  straight 
line  is  the  most  obvious  approximation  to  the  current -pressure  dependence 
on  the  logarithmic  plot. 
C.   Waves  of  Stratification  Excited  by  Acoustic  Waves 

Table  A-l  shows  the  general  increase  in  both  velocity  of  propagation 
of  the  waves  (up,  ur)  and  speed  of  the  individual  striations  (vD,  v  ) 
as  pressure  is  decreased.   Subscripts  p  and  r  are  commonly  used  to  desig- 
nate the  slow  and  fast  waves.   The  fast  wave  recorded  at  0.82  torr  may 
have  been  an  s  wave.  Wave  propagation  velocities  are  accurate  only  to 
16%  to  20%  because  amplification  is  large  in  many  cases,  and  at  times 
two  waves  have  probably  merged.   These  factors  make  it  difficult  to 
determine  the  maximum  of  the  wave  group  amplitude,  and  the  packets  are 
wide  enpugh  that  the  method  given  by  Krejci  and  Pekarek  in  I  34  cannot 
be  successfully  applied.   In  cases  where  the  sound  wave  was  propagated  in 
the  cathode  to  anode  direction,  this  acted  to  increase  the  apparent  group 
velocity,  and  the  acoustic  wave  also  gave  rise  to  the  striations.   Stria- 
tion velocity  measurements  are  accurate  to  5  to  10%,  and  wavelength 
accuracy  is  about  10%. 

Four  general  variations  in  parameters  of  the  wave  group  can  be  obtained 
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Table  A-l 

Characteristics  of  Acoustically   Excited 
Waves  of   Stratification 


Pressure 

i 

u 
P 

VP 

\ 

u 
r 

V 

r 

\ 

(torr) 

(mA) 

(m/s) 

(tn/s) 

(cm) 

(m/s) 

(m/s) 

(cm) 

4.0 

0.3 

36 

8.2 

2.4 

160 

35 

4.4 

4.0 

0.2 

X 

ft. 8 

2.0 

140 

28 

4.2 

3.3 

0.6 

64 

6.9 

2.3 

X 

33 

4.5 

2.7 

0.75 

50 

7.0 

2.4 

220 

39 

4.3 

2.7 

0.50 

65 

4.9 

2.1 

235 

28 

3.9 

2.5 

0.60 

49 

5.5 

2.1 

X 

32 

4.1 

2.2 

3.1 

67 

20/17 

2.8/2.7 

X 

X 

X 

2.2 

1.6 

57 

10.5 

2.6 

X 

31 

5.1 

2.2 

1.3 

50 

9.3 

2.5 

X 

X 

X 

2.2 

0.7 

61 

5.6 

2.3 

X 

X 

X 

1.9 

1.7 

61 

11 

2.6 

X 

X 

X 

1.6 

1.7 

82 

8.4 

2.6 

X 

X 

X 

1.4 

8.0 

133 

28 

2.9 

X 

X 

X 

1.4 

5.0 

108 

22 

2.9 

X 

X 

X 

1.4 

2.6 

100 

13 

2.7 

X 

X 

X 

1.4 

1.5 

105 

8.5 

2.5 

X 

X 

X 

1.3 

2.6 

84 

11 

2.6 

X 

X 

X 

1.0 

9.6 

140 

15 

3.0 

X 

X 

X 

0.93 

4.8 

X 

14 

2.6 

X 

X 

X 

0.82 

33 

220 

70 

3.0 

5000 

250 

5.4 

0.82 

10.9 

160 

36/24 

3.0/2.6 

X 

X 

X 
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Table 

A-l 

(continued) 

0.82 

10.9 

140 

24/36 

2.7/3.0 

0.77 

15 

140 

24 

2.9 

0.77 

10 

220 

30/21 

3.1/2.8 

0.77 

8.0 

220 

29/31 

3.0/2.7 

0.77 

5.0 

220 

21 

3.1 

X  X 

X  X 

X  X 
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from  Table  A-l:  (1)  wave  propagation  velocity  (u)  decreases  with  in- 
creasing pressure,  (2)  striation  velocity  (v)  increases  with  increasing 
current,  (3)  striation  wavelength  decreases  with  increasing  pressure, 
(4)  striation  wavelength  increases  with  increasing  current.  There  is  no 
theory  which  gives  satisfactory  values  of  u  and  v  as  functions  of  current 
and  pressure  for  conditions  other  than  the  zero  Debye  length  limit.   Early 
theoretical  works  of  Pekarek  related  these  velocities  to  plasma  para- 
meters through  the  relaxation  time    : 

2    u  +  v 

In  order  to  compare  these  predictions  of  relaxation  time  with  experi- 
mental observations,  acknowledge  of  plasma  parameters  over  the  experi- 
mental range  is  necessary.   The  measurements  of  plasma  parameters  which 
were  carried  out  do  not  cover  a  sufficient  range  of  current  and  pressure 
to  allow  this  comparison  to  be  made. 
D.   Irradiation  of  the  Discharge 
1.   Experimental  Arrangement 

Tube  I  was  filled  with  2.6  torr  neon,  sealed,  and  mounted 
parallel  to  a  longer  2.5  cm  diameter  tube,  with  a  centerline  separation 
distance  of  2.1  cm.   A  strip  of  black  paper  approximately  0.6  cm  wide 
was  attached  to  the  photomult iplier  housing.   This  encircled  tube  I  to 
prevent  light  of  the  illuminating  discharge  from  entering  the  slit.   A 
point  10  cm  from  the  anode  was  chosen  for  observing  the  effects  of  the 
irradiation  upon  the  light  intensity  at  different  values  of  discharge 
current.   Time-space  photographs  were  taken  with  27  cm  of  photomultiplier 
travel.   Current  in  the  illuminating  discharge  was  500  mA.   A  black 
cloth  at  times  covered  the  cathode  region  of  the  discharge  in  order  to 


100 


insure  that  all  effects  observed  were  the  result  of  changes  in  the  posi- 
tive column,  and  in  each  case  the  cloth  had  very  little  noticeable  effect 
on  the  changes  brought  about  by  irradiation. 

Peak  value  of  the  light  emitted  by  the  striations  in  the  illuminating 
discharge  was  found  to  increase  as  pressure  decreased.   Light  associated 
with  transitions  ending  on  metastable  states  is  the  important  quantity 
in  this  experiment,  not  total  light  emitted.   Eure  has  found  that  inten- 
sity of  the  8115  A  line  in  argon  increases  with  increasing  pressure.  I  18  J 

The  transition  involved  in  the  emission  of  this  line  corresponds  to  that 

o 
involved  in  the  strong  6402  A  line  of  neon,  with  the  principal  quantum 

number  being  four  in  argon,  three  in  neon.   Hence  a  similar  pressure  de- 
pendence would  be  expected,  and  it  was  found  that  illuminating  discharges 
had  a  larger  effect  upon  the  discharge  in  tube  I  as  the  pressure  was  in- 
creased.  However,  Druyvestyn  and  Penning  quote  data  which  show  that  the 
fraction  of  power  input  to  the  discharge  which  is  radiated  in  metastable- 
ended  transitions  increases  as  pressure  is  increased  from  0.1  to  1.0 
torr,  and  decreases  as  pressure  is  further  increased.   16   Kagan, 
Lyagushchenko  and  Khakhaev,  working  in  the  pressure  range  1-30  torr, 
also  found  a  monotonic  decrease  in  the  total  number  of  transitions  which 
end  on  metastable  states.   31 

The  increase  of  total  electrical  signal  from  the  multiplier  seen  in 
figure  >4—  6  may  be  attributed  to  the  change  of  discharge  color  from  red 
to  orange  as  the  pressure  is  reduced.   Since  the  S-4  cathode  of  the  1P21 
has  a  quantum  efficiency  which  is  rising  rapidly  in  the  color  interval 
red  to  yellow,  the  same  total  number  of  transitions  will  cause  the 
electrical  signal  to  be  larger  when  the  discharge  is  more  orange  (lower 
pressure) . 
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2.      Irradiation  Effects  on   Steady-State  Conditions 

Figure  4-6  shows   variation  of    light    intensity   and   potential 
across   the  discharge  tube   as   a  function  of  discharge  current.       (A   line 
drawing  is  given   in  figure  4-7  to   assist    in  explaining  this    figure.) 
In  part    (a)    of   figure   4-6,    the  discharge   is   not    illuminated.      Low  fre- 
quency   (slow)    striations    appear   as   current    is    increased   above   two  milli- 
amperes ,    and   high  frequency    (fast)   striations   appear  below  0.5   milli- 
amperes.      Potential   across   the  discharge  tube   increases   as   current  de- 
creases  until   the  fast    striations   appear  at  0.5  mA.      At  this   point   the 
potential    increases    less   rapidly,    indicating  that    less   power   is   dissipa- 
ted  in  the  discharge  when  the  fast    striations    are   present.      This    is    in 
accord  with  Pekarek's   finding  that   the  striated   state  of  the  column   is 
ire  conservative  state.      42       A   slight   oscillation   of  the  potential 


a  moi 


across  the  discharge  tube  in  the  presence  of  slow  striations  is  also 
evident. 

In  part  (b)  of  figure   4-6,  the  discharge  is  illuminated  by  a  500 
mA,  11.4  torr  discharge.   The  light  intensity  trace  shows  that  both 
striation  boundary  currents  have  been  increased;  the  lower  current  bound- 
ary from  0.5  to  0.85  mA,  and  the  upper  current  boundary  from  2.05  to 
2.15  mA.   Striations  appearing  at  the  lower  boundary  have  larger  ampli- 
tudes than  in  the  unilluminated  case,  while  striations  appearing  at  the 
upper  boundary  have  smaller  amplitudes.   The  striations  above  the  upper 
boundary  appear  to  have  maximum  amplitudes  which  are  less  well  defined 
than  those  in  part  (a)  of  this  photograph.   Moving  striations  were  pre- 
sent in  the  illuminating  discharge,  and  these  were  unstable  due  to  anode 
spot  oscillations.   These  oscillations  of  the  illuminating  discharge  may 
have  caused  variation  of  the  amplitude  of  the  slow  striations. 
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Figure  4~& .   Light  intensity  (increasing  to  right)  and 

voltage  drop  across  tube  I  at  2.6  torr.   Horizontal  scale: 
0.5  mA/div.   Vertical  scale:  20v/div.,  with  the  zero  of 
potential  offset.   Upper  -  no  irradiation.   Lower  -  500  mA 
in  irradiation  discharge. 
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Figure  4-7.   Light  intensity  and  potential  across  the  discharge 
tube  as  a  function  of  current.   This  diagram  is  a  key  to 
figure  4-14. 
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Illumination  of  the  discharge  by  light  from  a  discharge  in  the  same 
gas  reduces  the  population  of  metastable  states  by  promoting  atoms  in  the 
metastable  states  to  higher  energy  states  from  which  it  is  possible  to 
make  transitions  back  to  the  original  state  or  to  another  state.   Since 
the  population  of  the  metastable  states  is  reduced,  the  stepwise  ioniza- 
tion rate  is  reduced,  requiring  an  increase  in  the  rate  of  direct  ion- 
ization.  This  must  be  accomplished  through  an  increase  in  electron 
temperature,  which  in  turn  is  a  function  of  the  electric  field.   The 
trace  of  tube  potential  in  figure  4-6  (b)  clearly  shows  this  increase 
of  electric  potential  in  the  striation-free  and  slow  striation  regions. 
An  increase  is  also  present  in  the  current  range  of  the  fast  striations 
above  a  current  of  0.5  mA.   In  this  range  of  current  the  electric  poten- 
tial across  the  discharge  tube  not  only  increases  less  rapidly  with  de- 
creasing current,  but  is  found  to  decrease  immediately  below  the  stria- 
tion boundary  before  again  increasing  as  current  is  further  decreased. 

Working  with  this  concept  of  the  effect  of  the  illuminating  dis- 
charge, the  decrease  of  amplitude  of  the  slow  striations  and  increase  of 
the  upper  current  boundary  indicate  that  the  slow  striations  may  be 
associated  with  metastable  atoms.   Reduced  amplitude  of  these  striations 
is  explained  on  the  basis  of  the  lower  population  of  the  metastable 
levels.   Similarly,  the  increase  of  fast  striation  amplitude  accompany- 
ing the  increase  of  direct  ionization  under  illumination  indicates  that 
these  striations  are  associated  with  direct  ionization.   Pekarek  has 
also  stated  both  of  these  possibilities  from  his  relaxation  time  measure- 
ments and  from  his  irradiation  experiments.   47,  48 

Total  electrical  signal  from  the  photomultiplier  is  higher  when  the 
discharge  is  illuminated.   A  change  of  color  from  red  to  orange  is  also 
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observed,  and  the  higher  quantum  efficiency  of  the  photomultiplier  at 
shorter  wavelengths  (discussed  earlier)  may  account  for  this  increase. 

3.   Irradiation  Effects  on  Waves  of  Stratification 

Figure  *J—  8(a)  shows  a  slow  wave  of  stratification  in  tube  I 
for  a  discharge  current  of  2.0  mA.   Examination  of  the  striation  paths 
from  the  right  hand  side  of  the  page  (parallel  to  the  paths)  reveals 
convex  dispersion  in  the  cathode  half  of  the  column,  and  concave  disper- 
sion in  the  anode  half  of  the  column.   Part  (b)  of  the  figure  shows  the 
slow  wave  propagating  in  a  discharge  with  the  same  current.   In  this  case 
dispersion  remains  convex  throughout  the  column.   Amplification  of  the 
wave  is  also  decreased  by  illumination. 

Discussion  of  dispersion  showed  that  convex  dispersion  is  predicted 
by  zero  Debye-length  theories,  and  concave  dispersion  is  presently 
thought  to  be  due  to  a  quadratic  ionization  term  which  shifts  the  ion- 
ization maximum  in  the  direction  of  the  anode.   It  is  therefore  attribu- 
ted to  a  non-zero  Debye  length.   Since  the  Debye  length  exhibits  the 

2 
dependence  3  °<  T  /N  ,  an  increase  in  Te  and  a  decrease  in  Ne  both  act 

to  increase  ^D.  The  change  of  dispersion  observed  in  the  photograph  of 
figure  4-8  would  indicate  a  decrease  in  Debye  length  if  this  explanation 
were  correct. 

Assuming  that  the  irradiation  effects  can  be  accounted  for  on  the 
basis  of  changes  in  electron  density  and  temperature  alone,  the  follow- 
ing line  of  reasoning  applies: 

a)  E  is  increased 

b)  Te  is  increased 

c)  Since  j  =  <^E,  where 

j  =  current  density 
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Figure  *}■&.      Waves  of  stratification  in  tube  I  at  2.6  torr 
and  2.0  ma.   Time:   0.4  msec/div. ,  deflection:  27  cm. 
Upper  -  no  irradiation,  lower  -  500  mA  in  illumination 
discharge. 
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<f  -   plasma  conductivity 
and  &*  °*   N0. 
A  constant  j  requires  a  decrease  in  N  ,  therefore 
d)  Te/N0  is  increased,  and  A-^  is  increased. 
This,  however,  is  contrary  to  the  expected  result. 

Before  drawing  conclusions  as  to  the  validity  of  theoretical  ideas 
relating  to  dispersion,  it  must  be  pointed  out  that  insufficient  data 
is  available.  The  irradiation  effects  must  be  examined  for  several  cur- 
rents in  the  striation-free  region  to  insure  that  the  phenomenon  recorded 
is  a  general  occurrence.   Data  collected  during  this  work  indicate  that 
concave  dispersion  is  observed  less  often  than  convex  dispersion.  The 
problem  requires  further  study  before  a  statement  can  be  made. 
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APPENDIX  III 
OTHER  WORKS  ON  PLASMA-ACOUSTIC  WAVE  INTERACTIONS 

Suits  evidently  published  the  first  paper  related  to  this  subject 
in  1935.   81   Using  an  arc  discharge  as  his  main  discharge,  he  excited 
sound  waves  by  a  high-voltage  pulsed  spark  discharge  at  one  end  of  the 
arc  column,  and  detected  the  acoustic  waves  by  observing  the  variation 
in  electric  current  drawn  in  a  second  spark  discharge  placed  near  the 
other  end  of  the  column.   From  the  time  required  for  propagation,  he 
deduced  the  propagation  velocity  and  hence  the  neutral  gas  temperature 
in  the  arc.   Blum,  Dyer  and  Ingard  used  this  effect  of  gas  pressure  on 
current  in  their  construction  of  a  corona  discharge  for  use  as  an  acoustic 
transducer.  I  6  I  The  corona  is  a  special  case  of  a  glow  discharge,  and 
current  drawn  for  a  constant  voltage  is  a  function  of  mobility,  which 
in  turn  is  a  function  of  density.   They  indicated  that  this  was  a 
sensitive  microphone,  but  as  nothing  has  appeared  in  the  literature 
other  than  their  abstracts  of  papers  for  meetings  of  the  Physical  and 
Acoustical  Societies,  no  quantitative  information  is  available. 

Wojaczek  inaugurated  the  present-day  form  of  this  work  by  reporting 
modulation  of  the  light  intensity  of  the  positive  column  as  a  result  of 
a  high-voltage,  microsecond  pulse.   87 I   This  was  discovered  during 
investigations  of  the  positive  column  with  regard  to  waves  of  stratifi- 
cation.  He  found  that  the  effect  was  larger  when  propagation  was  in  the 
anode  to  cathode  direction.   An  estimate  of  the  damping  showed  it  to  be 
an  order  of  magnitude  larger  than  that  predicted  by  acoustic  wave  theory. 
Hayess  performed  this  experiment  for  the  express  purpose  of  examining 
the  damping  coefficient  as  a  function  of  current  and  pressure.   25  I 
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Working  in  xenon,   he  found  that   varying  the  current   had   little   effect 
on  the  damping   coefficient,    and   an   increase   in  pressure  resulted    in  a 
decrease   in  the  damping   coefficient.      Wall  damping  was   calculated  to  be 
the  most   important    influence   in  the  narrow   (1.55   cm  diameter)    tubes   used 
in  the  experiment.      This  was   an  order  of  magnitude   larger  than  damping 
due  to  collisions   and  thermal  conductivity.      Wojaczek  had   considered 
only  the  latter  two  mechanisms   in  his   calculations,    and  this   explains 
the  discrepancy  noted   above. 

Berlande,   Goldan   and  Goldstein  found  modulation  of  the   light  emit- 
ted by  an   afterglow  discharge   in  a  cryogenic   plasma.    J4J     They  attribute 
the   acoustic  wave  production  to  decay  of  electron  temperature,   but 
present   no   analysis   to   support   this   contention. 

Goldstein,   Roux,    and  Dayton   introduced   a   sound  wave    into  a  dis- 
charge tube  by  microwave  breakdown   and  heating  of  the  gas   at   one  end. 
I 221      This    sound  wave   was   detected  by  monitoring  the  variations  of 
quantities    in  three   spearate  devices:      (1)    light    intensity   in  an  after- 
glow plasma   located   at   some  distance   from  this    source, (2)   microwave 
power  transmitted   in  a  waveguide   in  which  the    same  afterglow  plasma  was 
used   as   a  dielectric   post, (3)    current    in  a  gas-filled  thermionic  diode. 
The   latter   is   explained    in  Chapter   II. 

Gentle   and  Ingard  used  a  high-voltage  pulsed   breakdown  for   sound 
wave  production,    and  biased  Langmuir  probes    into   the   ion  collection 
region  to  measure  the  sonic  velocity.      In  this  manner  they  determined 
gas   temperature  in  a  dc  discharge  maintained   in  a  flowing  gas.    I  2 X I 
Nygaard   also  used  a  high-voltage   pulse  breakdown  for  sound   wave  produc- 
tion,  observing  the  arrival  of  the   sound  wave  by  means  of   a  time  dis- 
placement  of  the  density-sensitive   Tonks-Dattner  resonances    in  an  after- 


110 


glow  plasma.    J  39  I 

Ingard   and  Gentle  have   shown  that   density  variations    in  the   neutral 
gas   will  cause  corresponding  variations    in  the   ion  and   electron  densi- 
ties.     27        Saxton's   work  confirmed  this.    172        Ingard's   later  work 
also   shows   that   acoustic  waves  may  be   amplified   in  the  discharge.    J  26  J 

Sodha  and  Palumbo  have  theoretically   analyzed    interaction  of 
acoustic  and  electromagnetic  waves    in  a  plasma.    I 75  j      They  note  that 
electromagnetic  waves    in  the  plasma  can  be  modulated   by  acoustic   waves 
in   such  cases.      Modulation  can  also  occur   upon  reflection  from  a  plasma- 
free-space  boundary  which   is   affected  by  an  acoustic   signal.      This   can 
become   important    in  upper  atmosphere  detection  of   acoustic  waves   gen- 
erated by  explosions   and  rockets. 

Strickler   and   Stewart   were   in  the  process  of   investigating  the 
effect   of   square-wave  modulation  of  the  discharge   in  connection  with 
moving  striations,   when  they  discovered  that   they  had   caused  the   column 
to  assume  a  permanent   kink.      771      This   was   attributed  to  the  excitation 
of  a  normal    longitudinal   acoustic  mode  of   the   column,    and   is  the  most 
striking  effect  observed   so   far. 
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